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ABSTRACT 

Simulations of the formation of large-scale structure predict that dark matter, low density highly 
ionized gas, and galaxies form 10-40 Mpc scale filaments. These structures are easily recognized in 
the distribution of galaxies. Here we use Lya absorption lines to study the gas in 30x6 Mpc filament 
at cz~3500 kms -1 , defined using a new catalogue of nearby (cx<10,000 kms -1 ) galaxies, which is 
complete down to a luminosity of about 0.05 L* for the region of space analyzed here. With HST 
spectra of 24 AGN we sample the gas in this filament. All of our sightlines pass outside the virial 
radius of any known filament galaxy. Within 500 kpc of the filament axis the detection rate is ~80%, 
but no detections are seen more than 2.1 Mpc from the filament axis. The width of the Lya lines 
correlates with filament impact parameter and the four BLAs in our sample all occur within 400 kpc of 
the filament axis, indicating increased temperature and/or turbulence. Comparing to simulations, we 
find that the recent Haardt & Madau (2012) extragalactic ionizing background predicts a factor 3-5 
too few ionizing photons. Using a more intense radiation field matches the hydrogen density profile 
within 2.1 Mpc of the filament axis, but the simulations still overpredict the detection rate between 
2.1 and 5 Mpc from the axis. The baryonic mass inside filament galaxies is 1.4xlO 13 M 0 , while the 
mass of filament gas outside galaxy halos is found to be 5.2 x 10 13 M 0 . 

Subject headings: galaxies: halos; intergalactic medium; quasars: absorption lines; ultraviolet: general 


1. INTRODUCTION 

The current paradigm for the formation of large-scale 
structure holds that after the Big Bang dark matter is 
concentrated by gravity into sheets, filaments and halos. 
Observational (e.g. Fukugita & Peebles 2006, Shull et 
al. 2012, Danforth et al. 2014) and theoretical (e.g. Cen 
& Ostriker 1999; Dave et al. 2001; Smith et al. 2011, 
Cen 2013) studies suggest that the baryons are carried 
along, with a small fraction (~10% at the present time) 
forming galaxies inside dark matter halos. The remain¬ 
ing baryons stay in the circumgalactic (CGM) and in¬ 
tergalactic (IGM) medium, with at z=0 about 30±10% 
in the form of photoionized H I at 10 4 K (Penton et al. 
2002; Lehner et al. 2007; Danforth & Shull 2008; Shull et 
al. 2012), and 40% to 60% at higher temperatures in the 
Warm-Hot Intergalactic Medium (WHIM). All these dif¬ 
ferent processes result in the Cosmic Web of dark matter, 
galaxy and gas filaments. 

The galaxies represent a small fraction of the baryons. 
As they are luminous and energetic, our understanding 
is relatively advanced. However, to fully understand the 
development of structure we should also understand the 
10 times more numerous (but more difficult to observe) 
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baryons outside of galaxies still in the IGM. 

Hydrodynamical simulations can be used to describe 
the evolution of Lya lines from high redshifts until the 
present. They predict the gradual disappearance of the 
Lya forest from z =3 to z= 0, the increase in temperature 
of the IGM over time, and the nature of the association 
of intergalactic gas with galaxies (Dave et al. 1999). Sim¬ 
ulations also predict that the optical depth of H I absorp¬ 
tion is closely tied to the underlying overdensity of dark 
matter (e.g. Croft et al. 1998; Schaye 2001), and give 
a specific expectation for the hydrogen column density 
profile perpendicular to the dark matter filaments. 

Observationally, the decrease in the number counts 
and column density of Lya absorbers over time pro¬ 
vides evidence for the theoretical picture of the IGM 
(e.g. Weymann et al. 1998; Kim et al. 2007, 2013; Rudie 
et al. 2013; Rahmati et al. 2013). Other statistical evi¬ 
dence comes from the spatial association between strong 
(log 7V(H l)>14) Lya absorbers and galaxies (Morris 
& Januzzi 2006; Ryan-Weber 2006; Stone et al. 2010) 
and from the power spectrum of the Lya absorbers (see 
e.g. Paschos et al. 2009; Lee et al. 2015 and references 
therein). 

To turn the simulated hydrogen density into a pre¬ 
diction for the observable H I column density requires 
a prescription for the intensity of the ionizing flux in 
the extragalactic background (EGB). Until recently, the 
most widely used prescription was provided by Haardt 
& Madau (2001), based on the spectra and number den¬ 
sity of quasars and galaxies and the evolution of those 
quantities over time. They updated their model in 2012 
(Haardt & Madau 2012), but this later version has much 
lower (a factor 3.7) ionizing flux at z=0. 

At z =0 Dave et al. (2010) compared the simulated col- 
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umn density distribution of Lya absorbers to the obser¬ 
vations of Lehner et al. (2007), showing a relatively good 
fit. However, with the updated version for the EGB given 
by Haardt & Madau (2012), and a larger sample of low 
redshift Lya absorbers (Danforth et al. 2014), Kollnreier 
et al. (2014) found a factor of five times too few photons 
are produced to ionize the Lya forest to the observed 
levels, by comparing the expected and observed column 
density distribution of Lya absorbers. In this paper we 
will confirm this deficit using a different method. The 
discrepancy has been further analyzed, with Khaire & 
Srianand (2015) arguing that the quasar contribution is 
higher by a factor of two, and additionally suggesting a 
4% escape fraction of Lyman continuum photons from 
star-forming galaxies to match the opacity inferred from 
the Lya forest. Shull et al. (2015) use different simula¬ 
tions to reduce the H I column densities, probably due 
to extremely strong feedback and argue that the Haardt 
& Madau (2012) quasar emissivity combined with ~5% 
escape fraction will produce the observed Lya forest. 

Most of the effort to understand absorption lines from 
the IGM has gone into work on the circumgalactic 
medium of galaxies (the gas within ~300 kpc), using en¬ 
semble studies of individual galaxy halos (see e.g. Morris 
et al. 1993; Lanzetta et al. 1995; Procliaska et al. 2006; 
Stocke et al. 2006; Wakker & Savage 2009; Rudie et al. 
2012; Tumlinson et al. 2013; Stocke et al. 2013 and many 
references in those papers). Yet, the most visually strik¬ 
ing prediction of the hydrodynamic simulations is the 
presence of large-scale filamentary structure. Such fila¬ 
mentary structure is clearly evident in the distribution 
of galaxies, but this gas has not been directly studied. 
Crosscuts through simulations (see e.g. Fig. 2 in Lukic 
et al. 2015) clearly show these filaments in the H I dis¬ 
tribution. 

The possible presence of filaments can be suggested by 
looking at the distribution of Lya absorbers in a single 
sightline. Morris et al. (1993) analyzed a deep galaxy 
sample within a degree of the 3C 273 sightline, combined 
with a shallower sample out to 10 Mpc. They found 
that the distribution of galaxies and absorbers as func¬ 
tion of redshift shows definite peaks (indicating the pres¬ 
ence of galaxy filaments), and that absorbers are not 
randomly distributed with respect to the galaxies; they 
also found absorbers that seemed to have no associated 
galaxies within 1 Mpc. In the same manner, Tripp et 
al. (1998) also showed that both absorbers and galaxies 
cluster around certain redshifts for the sightlines toward 
H 1821+643 and PG 1116+215. 

Studies of the transverse distribution of Lya absorbers 
with respect to the sightline are rare. Tejos et al. (2014) 
looked at the correlation between galaxies and absorbers 
in six fields toward QSOs with 2^0.7. However, they 
could only study the correlation function out to a few 
Mpc from each sightline. From this they concluded that 
50% of weak H I absorbers reside in what they call 
“galaxy voids”, i.e. no substantial dark matter halo is 
present within a few Mpc of the absorber. However, they 
did not really map the distribution of absorbers relative 
to the structures seen in the galaxy distribution. 

At redshifts beyond about 10,000 kms -1 mapping out 
the galaxies over large regions around individual sight¬ 
lines is observationally expensive. The galaxies are rel¬ 
atively faint but still need to be mapped over many de¬ 


grees, whereas instrumental limitations typically results 
in galaxy maps covering only about one degree around 
a sightline and larger-scale mapping is eschewed in favor 
of mapping multiple sightlines. At the lowest redshifts 
(c2<5,000 kms -1 or so), however, it is possible to get a 
mostly complete galaxy sample down to low luminosities 
(0.1 L* or lower) across most of the sky, since this lumi¬ 
nosity corresponds to a galaxy with r?z~ 17. Most of the 
sky has been covered this deeply, especially in regions 
covered by the footprint of surveys such as the Sloan 
Digital Sky Survey (SDSS) and the 6dF Galaxy Survey 
(6dFGS). 

The potential of this approach was first shown in a 
study by Narayanan et al. (2010). They combined Lya 
and O VI absorption at c,z~3000 kms -1 in the direction 
of the Seyfert galaxy Mrk290 to show the presence of 
gas at a temperature of 1.4xl0 5 K at an impact param¬ 
eter of 475 kpc to the nearest large galaxy, NGC5987. 
However, they also showed that this galaxy is located 
in a 30° (~30 Mpc) long T-shaped filament. With two 
other sightlines through the filament showing absorp¬ 
tion at similar velocities, while three off-filament sight¬ 
lines yielded non-detections, the interpretation of the ab¬ 
sorber toward Mrk 290 as originating in the outer halo of 
NGC 5987 was thrown into doubt. 

To follow up on this realization, we used the Cosmic 
Origins Spectrograph (COS) on the Hubble Space Tele¬ 
scope (HST) to obtain spectra for 17 additional AGN 
sampling the filament found by Narayanan et al. (2010). 
Supplementing these with seven archival sightlines, we 
describe in this paper how we used this data to develop 
a new approach to constraining the properties of the Lya 
forest that uses the spatial distribution of the absorbers 
in relation to the location of the dark matter/galaxy fil¬ 
aments. To do so we first need to set up several indepen¬ 
dent pieces: 

(a) Sect. 2 presents the galaxy filament found by 
Narayanan et al. (2010) and the details of the method 
that we used to objectively define the filament axis. 

(b) Sect. 3 presents the selection of targets and our defi¬ 
nition of a “filament impact parameter”. 

(c) In Sect. 4 we have a diversion in which we describe our 
method to correct the wavelength scale of COS spectra 

(d) Sect. 5 concentrates on the description of simulations 
from which H I column densities can be derived, needed 
to interpret our results. 

We combine all these items in Sect. 6, where we give the 
measurements on the spectra, discuss the galaxy/AGN 
impact parameters and the distribution of equivalent 
widths, column densities and line widths as function of 
separation to the filament axis, both in terms of just 
observations and in terms of comparisons to the sim¬ 
ulations. Specifically, we look at the influence of the 
assumed intensity of the extragalactic background radia¬ 
tion and compare the predictions for the run of detection 
fraction as function of filament impact parameter to the 
observations. 

2. THE GALAXY FILAMENT 
2.1. Galaxy sample and impact parameters 

In this section we describe the galaxy filament we study 
and discuss the method we use to derive a “filament im¬ 
pact parameter”, in addition to the usual “galaxy impact 
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Fig. 1.— Histogram of implied galaxy diameters, observed appar¬ 
ent magnitudes and derived luminosities and absolute magnitudes 
for the galaxies in the filament. This plot includes all 365 galaxies 
in the sky area shown in Fig. 2. Diameters are mostly based on the 
2MASS K-band angular diameters (see text), combined with the 
estimated galaxy distances. Luminosities follow from the empirical 
relation between luminosity and diameter found by Wakker & Sav¬ 
age (2009), which gives the luminosity to within a factor of 2. In 
the upper two panels vertical lines indicate galaxies with L=0,01, 
0.05, 0.5 and 1 L*. 

parameter” that is commonly used when studying inter- 
galactic absorption lines. 

We base our definition of the galaxy filament on the 
NASA Extragalactic Database (NED), first retrieving 
all information for all galaxies with cz<10,000 kms” 1 , 
finding their location on the sky, redshift, redshift- 
independent distance (if available), angular diameter, po¬ 
sition angle, inclination and type. This dataset will be 
described in more detail in a future paper (French et al., 
in preparation). 

Galaxy diameters were derived from angular diame¬ 
ters and individual galaxy distances. All angular di¬ 
ameter measurements were retrieved from NED. Diam¬ 
eters are preferentially taken from 2MASS K s “total” 
surface brightness extrapolation diameter measurements 
(for about 50% of all galaxies in our sample). For galax¬ 
ies with no 2MASS measurement, we plot different mag¬ 
nitudes against K s “total” values, and make a simple 
least squared linear fit between these magnitudes to es¬ 
timate K s . For ~20% of the sample SDSS (Sloan Digital 
Sky Survey) diameters are available, while for ~27% no 
diameter is given in NED. The remaining 3% of diame¬ 


ters are based on other surveys. To derive the 2MASS 
K s diameters, the “total” aperture radius rtot is de¬ 
fined to be the point at which the surface brightness 
extends down to 5 disk scale lengths (see Jarrett et al. 
2003 for a full description). The r to t value is derived as: 
rt.ot. = r' + a(lnM8) b , where r' is the starting point ra¬ 
dius (>5” - 10” beyond the nucleus, essentially beyond 
the PSF and nuclear or core influence) and a and b are 
the scale length parameters from a Sersic exponential 
function: / = f 0 exp(— r/a) (1 / b ). The fit extends to the 
point at which the mean surface brightness in the ellip¬ 
tical annulus has S/N<2. 

Figure 1 shows the distribution of diameters and ap¬ 
parent magnitudes of all galaxies in the region we ana¬ 
lyze (R.A. 208° to 276°, declination 41° to 72° and re¬ 
cession velocity 2900 to 4300 kms” 1 ; see Sect. 2.2). It 
also shows estimated luminosities and absolute magni¬ 
tudes, which are based on the empirical scaling between 
luminosity (L/Lfi) and diameter D that was found by 
Wakker & Savage (2009): logL = 2.31 logD—3.03. This 
relation yields luminosities to within a factor ~2. It gives 
L=0.05T* for D= 5.6 kpc and L=1L* for D=20.5 kpc. 
Using the fact that Mb=— 19.57 for a L* galaxy (Marzke 
et al. 1994), L=0.05L* corresponds to an apparent mag¬ 
nitude of 17.2 at a redshift cz=3500 kms. It would be 
preferable to use observed magnitudes and distances to 
derive galaxy luminosities, but the data in NED is ex¬ 
tremely inhomogeneous. We are working on a proper 
derivation of luminosity (French et al. in preparation), 
but preliminary numbers show that the distribution of 
luminosities will remain very similar. 

Fig. 1 shows that the number of galaxies increases 
down to a diameter of about 5 kpc (corresponding to 
a luminosity of 0.05~L„, i.e. about the luminosity of the 
SMC), with a break at about that dianreter/lunrinosity. 
We conclude that our galaxy sample is mostly complete 
above i=0.05L*. That this limiting luminosity is rela¬ 
tively low is due to the fact that this area of the sky is 
in the footprint of the Sloan Digital Sky Survey (SDSS); 
about half of the galaxies in the filament had their red- 
shifts first measured by this survey. 

We also estimate the virial radius of each galaxy, using 
the parametrization discussed by Stocke et al. (2013). 
This relates a galaxy’s luminosity to its virial radius. 
Fig. 1 of Stocke et al. (2013) shows a number of rela¬ 
tions between virial radius and luminosity, which they 
use to derive an a representative relation. This has a 
steeper slope at high luminosity than at low luminosity 
(L<0.1T*). Since almost all of our galaxies are brighter 
than 0.1 L* we use the plot in Stocke et al. (2013) to 
approximate their average relation between virial radius 
and luminosity as logi?. v j r = 0.3 log L + 2.25. Combin¬ 
ing this with the Wakker & Savage (2009) empirical re¬ 
lation between galaxy diameter and luminosity shows 
that a galaxy’s virial radius can be approximated as 
log l?vir=0.69 log D+ 1.24. 

2.2. Defining filament axes 

Figure 2 shows all galaxies with R.A. between 208° 
and 276°, declination between 41° and 72° and recession 
velocity between 2900 and 4300 kms” 1 . It is obvious that 
with this selection most galaxies (especially the ones with 
cz between ~3000 and ~3700 kms” 1 ) are distributed 
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Fig. 2. — Distribution of galaxies (colored circles) in the galaxy filament. The sizes of the galaxy circles are proportional to their area, 
while the color indicates their velocities, following the scale liar on the right. The labeled rectangles show the strip boxes used to calculate 
the filament axis. Box A is outlined with a thick line as a reminder that these boxes overlap by 50%. 


along two mostly-linear structures, which we interpret 
as a T-shaped galaxy filament. 

Several papers have been written about algorithms 
to automatically find filaments in galaxy datasets (e.g. 
Sousbie 2011, Cautun et al. 2013, Tempel et al. 2014, 
Chen et al. 2015, and references in these papers). How¬ 
ever, we do not have these codes available. Many are 
optimized for finding filaments in simulations, but not in 
real (messier) galaxy data, which have issues like varying 
luminosity limits and large distance uncertainties when 
using the Hubble constant to convert velocities to dis¬ 
tances. Until we can assess the usefulness of these various 
codes and algorithms, we proceed using the simpler, less 
general, but more intuitive method we developed our¬ 
selves, which is described below. 

1) First we draw two rectangular boxes on the sky map 
and divide these into 6 Mpc “strips” along the filament 
axis, using a distance of 50 Mpc, which is based on an 
estimated central velocity of 3500 kms -1 . These strips 
overlap by 50%, as shown in Fig. 2. 

2) Next, we make histograms of the velocities of the 
galaxies in each strip. These are shown on the left side of 
Fig. 3. This makes clear that there is a concentration of 
galaxies near ~3500 kms -1 , which is especially obvious 
in strips C through F and J through P. 

3) Using these histograms, we fit a gaussian to the ve¬ 
locity distribution of the galaxies in the filament, sepa¬ 
rately for each strip, determining a mean and a disper¬ 


sion. The mean velocities as function of location along 
the filament axis are then used to derive a linear rela¬ 
tion between the central filament velocity and position 
along the filament. The typical FWHM around the mean 
is about 600 kms -1 , which (combined with Fig. 3) we 
used to choose a velocity width of 1000 kms -1 for each 
strip. That is, we selected all galaxies in each strip within 
±500 kms -1 of the central filament velocity. The result¬ 
ing velocity ranges are indicated by the blue filled bins 
in Fig. 3. 

4) With the filament’s velocity structure now con¬ 
strained, we find the location of the filament axis by 
deriving the angular separation of each galaxy from the 
long edge of the guide box, i.e., the outside boundary of 
the strips shown in Fig. 2. This was done for each strip 
separately, as shown in the right-hand panels of Fig. 3. 
We then calculated the centroid and dispersion of this 
perpendicular distribution (using a gaussian fit if there 
are more than 12 galaxies, and a moment calculation if 
there are fewer). Combining the center of the strip along 
the long axis with the centroid of the galaxy distribution 
perpendicular to the axis defines a point on the filament 
axis. The axis itself is defined by connecting these points, 
creating 3 Mpc long axis segments. 

We note that the lower filament seen in Fig. 2 shows 
a slight velocity gradient, with the central filament ve¬ 
locity increasing from 3300 kms -1 in the lower right to 
3850 kms -1 in the upper left. The upper segment has 
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Fig. 3a. — The red histograms in the left panels give the ve¬ 
locity distribution of the galaxies in each strip (same label as in 
Fig. 2, given in the upper left corner of each panel). The filled blue 
bins show the galaxies considered part of the filament after ap¬ 
plying both the velocity and edge-separation selection, i.e. within 
±500 kms -1 of the central velocities and inside a selected edge- 
separation range. The red dashed vertical lines in the left panel 
give the velocity range included in Fig. 2, while the blue dashed 
vertical lines show the range of velocities used to select the galaxies 
in the filament. The right panels gives the distribution of distance 
to the edge of the strip for galaxies in the velocity range selected 
for each strip. The red histograms indicate galaxies in the proper 
velocity range, but not in the proper range of edge-separations. 
For strips A-H the edge of the strip runs along the bottom, from 
(RA,Dec)=(223,41) to (261,65), while for strips I-P this edge runs 
from (270,54) to (218,64) (see Fig. 2). The distribution of angular 
separations from this edge is determined by fitting by a gaussian 
if there are 12 or more galaxies, by the first and second moment if 
there are fewer than 12. The resulting centers are used to define 
the filament axis; these fits are shown by the solid blue line. 

a much smaller gradient, going from 3750 kms -1 on the 
left to 3550 kms -1 on the right. The best single velocity 
range that includes all galaxies in the filament runs from 
2900 to 4300 kms -1 , which is the range shown in Fig. 2. 
We note, however, that for our analysis we include the 
effects of the velocity gradient. 

Finally, we note that the filament is mostly parallel 
to the plane of the sky (which is why it can be easily 
recognized in a map selecting galaxies in a limited ve¬ 
locity range). The 550 kms -1 gradient corresponds to a 
change in distance of ~8 Mpc, contrasting with a pro¬ 
jected length of ~27 Mpc. 

3. OBSERVATIONS 
3.1. Datasets 

To sample the selected galaxy filament shown in Fig. 2, 
we constructed a set of 24 targets. First, we corre¬ 


Fig. 3b. — continued. 

lated the Veron-Veron QSO catalog with the GALEX 
database, to find all AGNs with UV flux larger than 
1.5xl0 -15 erg cm -2 s -1 A -1 . This gave a list of 75 ob¬ 
jects in the region of interest. We selected most of the 
brighter ones, but kept a few fainter targets located in 
strategic directions. The final set of 17 targets were ob¬ 
served using the Cosmic Origins Spectrograph (COS) on 
the Hubble Space Telescope (HST, program 12276) which 
has a resolution of about 20 kms -1 , allowing us to re¬ 
solve the Lya absorbers. That sample is supplemented 
by seven targets in the same region for which archival 
data are available, mostly also obtained using COS, but 
in one case (3C 351.0) with the Space Telescope Imag¬ 
ing Spectrograph (STIS). The STIS data are described 
by Tripp et al. (2008), and we use the same dataset for 
3C 351.0 as they did. The COS instrument is described 
in Green et al. (2012) and we use the xld spectra avail¬ 
able from the Multimission Archive at Space Telescope 
(MAST), with modifications (see Sect. 4). 

The locations and program information for each target 
are given in Table 1. The observations in program 12276 
were taken with central wavelengths set to 1291A and 
1327A, resulting in continuous wavelength coverage from 
~1135 to ~1465A. For six targets the exposure time is 
about 2 ks, while for ten targets two orbits were used 
(giving 5 ks exposures). One target was exposed for 3 
orbits (8.4 ks). When observing for multiple orbits, we 
also used two different FP-POS settings in the second or¬ 
bit (central wavelength 1327 A). The other seven targets 
were part of eight other programs, using a variety of ex¬ 
posure times and grating settings. The different central 
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wavelengths allow an assessment of detector glitches and 
fixed-pattern noise, as the same wavelength will fall on 
two or three different places on the detector. The CAL- 
COS pipeline (v2.19.7) was used to process the raw data, 
producing one-dimensional extracted spectra. 

3.2. AGN sampling the filament 

Figure 4 shows the directions to the 24 relevant AGN 
relative to the galaxy filament found in Sect. 2.2. This 
figure also shows cases where no Lya absorption is found 
in the velocity range 2900 to 4300 kms -1 as open stars 
and the 10 sightlines with detected Lya shown as colored 
stars. These measurements will be discussed in Sect. 6. 

Table 2 presents the properties of the galaxies near 
each AGN sightline. It includes all galaxies within the 
relevant velocity range (2400 to 4800 kms -1 , see below) 
that lie within 1 Mpc of the AGN sightline; if there are 
fewer than three such galaxies, it gives the three galaxies 
with the smallest impact parameter. 

For AGN toward which a Lya line is seen (see Sect. 6), 
Col. 2 shows the velocity of the absorber. Column 4 
then gives the difference between that velocity and the 
filament velocity. For example, 3C 351.0 has a Lya line 
at 3597 kms -1 , while the nearest axis segment (at im¬ 
pact parameter 569 kpc) is centered at 3600 kms -1 . This 
implies the absorption is offset by —3 kms -1 from the 
filament velocity. In the case of galaxies near the sight¬ 
line, this column also gives the difference between v(Lya) 
and v (galaxy), but only when that difference is less than 
±400 kms -1 and either the impact parameter is less than 
300 kpc or the ratio of impact parameter to virial radius 
(Col. 8) is less than two. That is, a non-blank entry in 
this column for a galaxy row means that in a conven¬ 
tional approach the Lya absorption could be associated 
with the galaxy. 

We note that a galaxy’s halo probably does not have a 
well-defined edge. Gas beyond the virial radius may be 
falling in if it happens to be moving radially. Other gas 
beyond the virial radius may be moving mostly trans¬ 
versely and never reach the galaxy. In order to better 
separate possible halo gas from possible filament gas, we 
use the arguments in the papers by Oort (1970), Mailer 
& Bullock (2004) and Shull (2014). They find that for 
a Milky Way sized galaxy (~L*) the boundary between 
gas falling in and being part of the IGM is at about 230 
kpc. In particular, Shull (2014) discusses the difference 
between this “accretion radius”, a “gravitational radius” 
and the virial radius. The approximation formula given 
in Sect. 2.1 gives 140 kpc as the virial radius for an L* 
galaxy. Thus a number on the order of 1.5 R„, r is justi¬ 
fiable as giving the border between gas in a galaxy halo 
and gas in the filament. 

Using this shows that there are no L>0.05 L* galaxies 
within these limits that can be associated with the Lya 
detections toward 4C 63.22, FBS 1526±659, H 1821±643, 
Mrk 486, RBS 1483, RBS 1503, RX J1500.5±5517 and 
RX J1717.6±6559, so conventionally these would be 
called “void absorbers”. 

There are galaxies within the 400 kms -1 , and 
300 kpc or 2 R t ,j r limit near the Lya absorbers seen 
toward 3C 351.0, Mrk 290, Mrk 876, RX J1608.3+6814, 
SBS 1537±577 and SBS 1551±572. However, only to¬ 
ward SBS 1537±577 is the nearest (small) galaxy within 


1.5 virial radii, while for the other sightlines the ratio of 
impact parameter to virial radius lies between 1.7 and 
2.3. I.e., these sightlines still pass rather far from the 
galaxies and they do not really sample the circumgalac- 
tic medium of these galaxies. 

Thus, using generous standard criteria 
(Aw<400 kms -1 , p<300 kpc) at most four of the 
thirteen Lya absorbers would be associated with a 
galaxy. Using even more lenient criteria (including 
galaxies up to 2 virial radii from the sightline) six of the 
detections could be associated. Using more physically 
plausible criteria (Au<400 kms -1 and / o/R vir <1.5) only 
one Lya absorber samples a galaxy halo. Therefore, 
basically all Lya absorbers in our sample are unlikely 
to be associated with the halos of galaxies brighter than 
the SMC (~0.05L*). 

We can estimate how likely it is that our sample of 
randomly placed sightlines passes inside the virial radius 
of one of the galaxies. Taking the implied virial radii 
of all galaxies in our filament and adding up the total 
sky area that is covered gives a value of 14 Mpc 2 . Since 
the lower branch of the filament covers about 18° x 7° 
(~15.5x6 Mpc), while the upper branch covers about 
15°x7° (^13x6 Mpc), the total filament area is about 
170 Mpc 2 . Thus, only about 8% of the filament area is 
covered by galaxy halos, and only 1 in 12 Lya absorbers 
should be associated with a galaxy. This is indeed what 
we find. 

As we will now argue, our Lya detections are instead 
likely to be associated with the IGM in the galaxy fila¬ 
ment. To quantify the relation between detections and 
the filament, we define a new quantity: the “filament 
impact parameter” (p(fil)) as the separation between an 
AGN sightline and the nearest axis segment. This is the 
product of the angular separation scaled by the distance 
corresponding to the segment’s recession velocity, i.e. the 
central velocity of each filament strip, found in the man¬ 
ner described above. The resulting impact parameters 
are given in Col. 4 Table 2 under the entry “Filament 
Axis” after each AGN’s name. 

4. COS DATA PREPARATION 
4.1. Correcting the COS wavelength scale 

As we were preparing COS spectra for a study of O VI 
absorbers at high S/N ratio (which turned into Savage et 
al. 2014) we discovered that in many cases spectral lines 
in different spectra of the same target did not align, with 
misalignments of up to ±40 kms -1 , and with the mis¬ 
alignment varying as function of wavelength. A summary 
of the COS wavelength calibration procedure is given by 
Oliveira et al. (2010). It is based on spectra of 6 bright 
targets taken in program 11474 and 11487, tying absorp¬ 
tion lines in their spectra to STIS-E140H data of the 
same targets. The workshop paper states that dispersion 
relations were derived only using the FP-POS=3 setting 
and that a linear dispersion relation was used for the 
G130M and G160M gratings. Thus, the offsets may be 
due to the fact that there is no explicit calibration spec¬ 
trum for every combination of central wavelength and 
FP-POS settings, so that small distortions in the im¬ 
age on the detector are not accounted for by CALCOS, 
which assumes the same dispersion relation for different 
FP-POS settings at the same central wavelength. Offsets 




Fig. 4. — Distribution of galaxies (colored circles) and AGN (open and colored stars) in the galaxy filament. The sizes of the galaxy 
circles are proportional to the square of their diameters (i.e. area), while the color indicates their velocities, following the scale bar on the 
right. AGN toward which Lya absorption is seen are shown by filled colored stars, while non-detections are indicated by open stars. The 
axes of the two filament segments are shown by the wiggly black lines. See the text for the derivation of these axes. 
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Fig. 5.— In each panel red points show the cross-correlation offsets for individual absorption lines between a spectrum and the chosen 
reference spectrum (blue points in panel 1 are for lines with high optical depth for which the centroid determination is more uncertain). The 
error bars indicate the uncertainties. The left panel gives the differences between the centroids of the ISM lines fitted in the 21-cm spectra 
vs. the combined UV spectra. The red lines are least-squares fits through these points. The blue lines in the spectrum-vs-spectrum panels 
show the final offsets implied for each spectrum, found by combining the relative offsets of the UV spectra with the offset needed to align 
the spectrum with the 21-cm data. The green line shows the offset needed to bring the nominally calibrated heliocentric velocities in line 
with LSR velocities. The label (e.g. Ibi603o51:1327:3(l)) contains, respectively, the observation ID, central wavelength, FP-POS setting 
and exposure ordinal number. Clearly, spectrum with 1327:3 show shifts relative to spectra with 1291:3 that vary by up to 60 kms -1 
across one side. 
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TABLE 1 
Observations 


Target 

(1) 

R.A. 

(2) 

Dec. 

(3) 

z 

(4) 

Program 

(5) 

PI 

(6) 

Obs. ID 
(7) 

Obs. date 
(8) 

Texp [fc«] 

(9) 

3C309.1 

14 59 07.6 

+71 40 20 

0.9050 

12486 

Bowen 

LBP240 

2012 03 27 

9.4 

3C309.1 



0.9050 

12486 

Bowen 

LBP245 

2012 03 29 

9.4 

3C351.0 

17 04 41.4 

+60 44 31 

0.3719 

8015 

Jenkins 

057901 

1999 06 27 

19.8 

3C351.0 



0.3719 

8015 

Jenkins 

057903 

1999 06 29 

14.4 

3C351.0 



0.3719 

8015 

Jenkins 

057902 

2000 02 10 

18.6 

3C351.0 



0.3719 

8015 

Jenkins 

057904 

2000 07 25 

26.5 

4C63.22 

15 23 45.9 

+63 39 24 

0.2040 

12276 

Wakker 

LBI601 

2011 08 14 

1.9 

FBS1526+659 

15 27 28.5 

+65 48 10 

0.3450 

12276 

Wakker 

LBI603 

2011 05 02 

2.0 

H1821+643 

18 21 57.2 

+64 20 36 

0.2970 

11484 

Hartig 

LABO06 

2009 07 26 

0.6 

H1821+643 



0.2970 

12038 

Green 

LBGL33 

2012 07 06 

11.5 

Kaz447 

17 03 28.9 

+61 41 09 

0.0773 

12276 

Wakker 

LBI618 

2011 10 05 

5.2 

Mrk290 

15 35 52.3 

+57 54 09 

0.0296 

11524 

Green 

LB4Q02 

2009 10 28 

3.9 

Mrk486 

15 36 38.3 

+54 33 33 

0.0389 

12276 

Wakker 

LBI607 

2011 12 18 

5.0 

Mrk817 

14 36 22.1 

+58 47 40 

0.0315 

11505 

Noll 

LACD01 

2009 08 04 

2.2 

Mrk817 



0.0315 

11524 

Green 

LB4Q01 

2009 12 28 

1.3 

Mrk876 

16 13 57.2 

+65 43 10 

0.1290 

11686 

Arav 

LB4F05 

2010 04 09 

3.1 

Mrk876 



0.1290 

11524 

Green 

LB4Q03 

2010 04 08 

9.5 

PG1626+554 

16 27 56.2 

+55 22 32 

0.1330 

12029 

Green 

LBGB01 

2011 06 15 

3.3 

RBS1483 

15 19 21.6 

+59 08 24 

0.0781 

12276 

Wakker 

LBI612 

2011 06 25 

1.9 

RBS1503 

15 29 07.5 

+56 16 06 

0.0990 

12276 

Wakker 

LBI614 

2010 11 26 

2.0 

RXJ1500.5+5517 

15 00 30.8 

+55 17 09 

0.4048 

12276 

Wakker 

LBI602 

2011 05 11 

8.4 

RX J1503.2+6810 

15 03 16.5 

+68 10 06 

0.1140 

12276 

Wakker 

LBI609 

2010 12 31 

1.9 

RXJ1508.8+6814 

15 08 52.8 

+68 14 07 

0.0586 

12276 

Wakker 

LBI608 

2010 12 05 

1.9 

RX J1608.3+6018 

16 08 20.5 

+60 18 28 

0.1780 

12276 

Wakker 

LBI610 

2010 11 21 

5.2 

RX J1717.5+6559 

17 17 38.0 

+65 59 39 

0.2927 

12276 

Wakker 

LBI619 

2011 06 22 

5.4 

SBS1458+535 

14 59 49.6 

+53 19 09 

0.3380 

12276 

Wakker 

LBI611 

2011 10 21 

5.0 

SBS1503+570 

15 04 55.6 

+56 49 20 

0.3589 

12276 

Wakker 

LBI617 

2011 10 19 

5.2 

SBS1521+598 

15 21 53.8 

+59 40 21 

0.2862 

12276 

Wakker 

LBI613 

2011 06 12 

5.1 

SBS1537+577 

15 38 10.0 

+57 36 13 

0.0734 

12276 

Wakker 

LBI606 

2011 10 19 

5.2 

SBS1551+572 

15 52 32.7 

+57 05 17 

0.3660 

12276 

Wakker 

LBI615 

2011 03 13 

5.1 

SBS1624+575 

16 25 26.5 

+57 27 27 

0.0670 

12276 

Wakker 

LBI616 

2011 10 17 

5.1 
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TABLE 2 

Galaxies between 2400 and 4800 kms -1 within 1 Mpc of each sightline 


Target 

Galaxy 

cz 

Av 

P 

Diameter 

L/L* 

R-vir 

p/ Rv ir 

Type 



[kms -1 ] 

km s —1 

[kpc] 

[kpc] 





(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 


3C309.1 

Filament Axis 

3600 


3908 






3C309.1 

IC1110 

3373 


3854 

24.6 

0.86 

158 

24.3 

Sa 

3C309.1 

UGC09734 

3341 


4756 

19.9 

0.55 

137 

34.7 

Sm 

3C309.1 

CGCG338-038 

3827 


4737 

15.9 

0.34 

117 

40.3 


3C351.0 

Lya @ 3597,3459 









3C351.0 

Filament Axis 

3600 

-3,-141 

569 






3C351.0 

Mrk0892 

3617 

-20,-158 

179 

12.9 

0.22 

101 

1.8 

Pair 

3C351.0 

NGC6307 

3057 


284 

22.7 

0.72 

149 

1.9 

PSBSOP. 

3C351.0 

NGC6306 

2973 


271 

19.3 

0.51 

134 

2.0 

.SBS2P* 

3C351.0 

NGC6310 

3419 


406 

34.5 

1.76 

200 

2.0 

.S..3* 

3C351.0 

NGC6292 

3411 


375 

22.4 

0.70 

148 

2.5 

.S..4.. 

3C351.0 

SDSSJ170327.95+610631.5 

3313 


333 

9.2 

0.11 

81 

4.1 


3C351.0 

2MASXJ17071270+6055144 

3099 

360 

280 

3.6 

0.01 

42 

6.5 

HII 

3C351.0 

SBS1700+603 

3736 


616 

9.6 

0.12 

83 

7.4 


3C351.0 

SDSSJ171138.94+604341.8 

3855 


828 

8.3 

0.08 

75 

11.0 


3C351.0 

SDSSJ171140.34+604115.6 

3350 


727 

5.5 

0.04 

57 

12.7 


3C351.0 

UGC10745 

3059 


700 

2.7 

0.01 

35 

19.8 

Sdm 

4C63.22 

Lya @ 2420 









4C63.22 

Filament Axis 

3600 

-1180 

2333 






4C63.22 

KHG1-C07 

2548 


2413 

1.0 

0.00 

3 

679.2 


4C63.22 

SBS1543+593 

2698 


3564 

1.0 

0.00 

3 

1003.3 

dwarf 

4C63.22 

KHG1-C09 

2848 


3593 

1.0 

0.00 

3 

1011.5 

Comp 

FBS1526+659 

Lya @ 3476 









FBS1526+659 

Filament Axis 

3600 

-124 

369 






FBS1526+659 

KHG1-C07 

2548 


3693 

1.0 

0.00 

3 

1039.8 


FBS1526+659 

KHG1-C09 

2848 


4593 

1.0 

0.00 

3 

1293.1 

Comp 

FBS1526+659 

SBS1543+593 

2698 


4821 

1.0 

0.00 

3 

1357.3 

dwarf 

H1821+643 

Lya @ 2825,4087 









H1821+643 

Filament Axis 

3600 

-775,487 

3249 






H1821+643 

NGC6636NED01 

4393 


2491 

40.2 

2.4 

222 

11.2 

,s?.... 

H1821+643 

NGC6687 

3374 


4239 

28.6 

1.18 

175 

24.1 

SAd 

H1821+643 

NGC6701 

3965 


4369 

29.6 

1.27 

180 

24.3 

PSBS1.. 

Kaz447 

Filament Axis 

3600 


184 






Kaz447 

NGC6310 

3419 


762 

34.5 

1.76 

200 

3.8 

.S..3* 

Kaz447 

NGC6292 

3411 


669 

22.4 

0.70 

148 

4.5 

.S..4.. 

Kaz447 

NGC6307 

3057 


826 

22.7 

0.72 

149 

5.5 

PSBSOP. 

Kaz447 

SDSSJ170327.95+610631.5 

3313 


485 

9.2 

0.11 

81 

6.0 


Kaz447 

NGC6306 

2973 


817 

19.3 

0.51 

134 

6.1 

.SBS2P* 

Kaz447 

Mrk0892 

3617 


949 

12.9 

0.22 

101 

9.3 

Pair 

Kaz447 

2MASXJ17071270+6055144 

3099 


699 

3.6 

0.01 

42 

16.3 

HII 

Kaz447 

UGC10745 

3059 


673 

2.7 

0.01 

35 

19.0 

Sdm 

Mrk290 

Lya; @ 3089 









Mrk290 

Filament Axis 

3600 

-511 

97 






Mrk290 

NGC5987 

3010 

79 

475 

46.4 

3.3 

245 

1.9 

.S..3.. 

Mrk290 

2MASXJ15351422+5730529 

3092 


318 

8.5 

0.09 

76 

4.1 

HII 

Mrk290 

SDSSJ153802.75+573018.3 

3525 


446 

10.4 

0.14 

87 

5.1 

Sd(f) 

Mrk290 

CGCG297-017 

3282 


506 

10.4 

0.14 

87 

5.8 


Mrk290 

SDSSJ153733.00+583447.8 

2932 


542 

10.7 

0.15 

89 

6.0 

Sc(f) 

Mrk290 

SBS1533+574A 

3348 


566 

6.8 

0.06 

66 

8.6 

HII 

Mrk290 

SDSSJ153040.88+575301.0 

2896 


519 

5.7 

0.04 

58 

8.9 


Mrk290 

SDSSJ153742.05+570506.4 

3469 


762 

7.2 

0.06 

68 

11.2 


Mrk290 

SBS1540+576 

3717 


760 

6.7 

0.05 

65 

11.7 


Mrk290 

SDSSJ153706.68+585651.5 

2989 


819 

6.6 

0.05 

64 

12.6 


Mrk290 

SBS1533+574B 

3429 


579 

3.9 

0.02 

44 

12.9 

HII 

Mrk290 

2MASXJ15335796+5650509 

3260 


915 

3.6 

0.01 

42 

21.5 


Mrk290 

SDSSJ152956.69+582635.9 

2908 


718 

2.5 

0.01 

33 

21.7 


Mrk486 

Lya <S) 4386 









Mrk486 

Filament Axis 

3600 

786 

1910 






Mrk486 

SBS1543+593 

2698 


3342 

1.0 

0.00 

3 

940.8 

dwarf 

Mrk486 

KHG1-C07 

2548 


3841 

1.0 

0.00 

3 

1081.3 


Mrk486 

KHG1-C09 

2848 


4615 

1.0 

0.00 

3 

1299.3 

Comp 

Mrk817 

Filament Axis 

3650 


6208 






Mrk817 

[DYC2005]443 

2998 


4679 

1.0 

0.00 

3 

1317.4 

Blue 

Mrk817 

RCS06100303673 

3298 


4717 

1.0 

0.00 

3 

1327.9 


Mrk817 

RCS06100304426 

3298 


4742 

1.0 

0.00 

3 

1335.1 
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TABLE 2 

Galaxies between 2400 and 4800 kms -1 within 1 Mpc of each sightline 


Target 

Galaxy 

cz 

Av 

P 

Diameter 

L/L, 

Rv'iy 

P/Rvir 

Type 



[kms -1 ] 

km s -1 

[kpc] 

[kpc] 





(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 


Mrk876 

Lya @ 3476 









Mrk876 

Filament Axis 

3600 

-124 

540 






Mrk876 

UGC10294 

3504 

-28 

274 

25.6 

0.93 

162 

1.7 

Sm? 

Mrk876 

UGC10376 

3246 


797 

19.3 

0.51 

134 

5.9 

Sm: 

Mrk876 

NGC6135 

3644 


685 

13.9 

0.25 

107 

6.4 


PG1626+554 

Filament Axis 

3550 


5684 






PG1626+554 

NGC6258 

3064 


4753 

22.4 

0.70 

148 

31.9 

.E. 

PG1626+554 

KHG1-C09 

2848 


4912 

1.0 

0.00 

3 

1383.0 

Comp 

PG1626+554 

SBS1543+593 

2698 


4919 

1.0 

0.00 

3 

1384.9 

dwarf 

RBS1483 

Ly a @ 2726 









RBS1483 

Filament Axis 

3600 

-874 

2260 






RBS1483 

NGC5894 

2466 


839 

33.6 

1.66 

196 

4.3 

.SB.8?. 

RBS1483 

SDSS J 151827.06+582658.6 

2555 


469 

7.8 

0.07 

72 

6.5 


RBS1483 

Mrk0847 

2540 


658 

10.2 

0.13 

86 

7.6 

,S?.... 

RBS1483 

SDSSJ152209.95+583819.5 

3349 


535 

5.8 

0.04 

59 

9.0 


RBS1483 

UGC09837 

2657 


880 

7.3 

0.06 

69 

12.7 

SAB(s)c 

RBS1483 

SDSS J151917.62+603435.9 

2535 


953 

4.8 

0.03 

51 

18.4 


RBS1503 

Lya @ 3269,3306 









RBS1503 

Filament Axis 

3600 

-331,-294 

306 






RBS1503 

NGC5965 

3412 


677 

55.6 

4.9 

277 

2.4 

.S..3.. 

RBS1503 

Mrk0482 

3357 


640 

13.2 

0.23 

103 

6.2 


RBS1503 

CGCG274-047 

3355 


517 

9.3 

0.11 

81 

6.3 


RBS1503 

2MASXJ15332453+5636315 

3553 


624 

10.2 

0.13 

86 

7.2 


RBS1503 

Mrk0481 

3298 


726 

8.9 

0.10 

78 

9.2 

E 

RBS1503 

MCG+09-26-001 

3366 


793 

9.1 

0.10 

80 

9.9 


RBS1503 

SDSSJ153325.63+564156.5 

3362 


636 

6.1 

0.04 

61 

10.4 


RBS1503 

SDSS J153007.14+553432.8 

3356 


612 

5.1 

0.03 

54 

11.3 


RBS1503 

SBS1524+554 

3409 


978 

8.5 

0.09 

76 

12.8 

BlueCG 

RBS1503 

SDSS J153558.83+564108.2 

3213 


859 

6.1 

0.04 

61 

14.0 

BLAGN 

RBS1503 

2MASXJ15335796+5650509 

3260 


744 

3.6 

0.01 

42 

17.5 


RXJ1500.5+5517 

Lya @ 3592 









RXJ1500.5+5517 

Filament Axis 

3650 

-58 

2031 






RXJ1500.5+5517 

UGC09737 

3367 


911 

16.1 

0.35 

118 

7.7 

S 

RXJ1500.5+5517 

SDSS J150621.06+550842.1 

3374 


738 

7.7 

0.07 

71 

10.3 


RXJ1500.5+5517 

SDSS J150804.21+551954.0 

3400 


948 

7.4 

0.07 

69 

13.6 


RXJ1500.5+5517 

SDSS J150654.40+553218.2 

3295 


805 

4.8 

0.03 

51 

15.5 


RXJ1500.5+5517 

SDSS J145718.28+543105.9 

3216 


750 

2.8 

0.01 

36 

20.7 


RXJ1503.2+6810 

Filament Axis 

3600 


1016 






RX J1503.2+6810 

IC1110 

3373 


1003 

24.6 

0.86 

158 

6.3 

Sa 

RXJ1503.2+6810 

UGC09855 

3480 


2538 

22.2 

0.69 

147 

17.2 

Im: 

RXJ1503.2+6810 

CGCG297-009 

2500 


4633 

1.0 

0.00 

3 

1304.4 

Sd(f) 

RXJ1508.8+6814 

Filament Axis 

3600 


1202 






RX J1508.8+6814 

IC1110 

3373 


799 

24.6 

0.86 

158 

5.0 

Sa 

RX J1508.8+6814 

UGC09855 

3480 


2249 

22.2 

0.69 

147 

15.2 

Im: 

RXJ1508.8+6814 

CGCG297-009 

2500 


4624 

1.0 

0.00 

3 

1301.9 

Sd(f) 

RXJ1608.3+6018 
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TABLE 2 

Galaxies between 2400 and 4800 kms -1 within 1 Mpc of each sightline 


Target 

Galaxy 

cz 

At; 

P 

Diameter 

L/L* 

Rvir 

p/ Rvir 

Type 



[kms -1 ] 

km s -1 

[kpc] 

[kpc] 





(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 
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can also occur if the target is not at precisely the same 
spot in the aperture during different visits; this can not 
be accounted for in the calibration code even in principle. 
Since the offsets can be as large as a resolution element, 
this produces an artificial smearing of the lines in the 
combined spectra. Moreover, absorption lines in a single 
absorption-line system can appear misaligned. 

Other authors have also noticed this issue. E.g., Sav¬ 
age et al. (2011), Meiring et al. (2011) and Tumlinson 
et al. (2011). They fix the offsets by cross correlating 
each spectral line separately, determining a centroid and 
adjusting the wavelengths for each line in each exposure 
individually. Different lines in the same absorption sys¬ 
tem are also explicitly forced to be aligned, although 
that is of course an assumption. Danforth et al. (2014) 
released a code in which they assume that there is a 
constant (wavelength-independent) offset between expo¬ 
sures, which is determined by a blind cross-correlation of 
a single 10 A wide spectral region around a strong inter¬ 
stellar line, separately for each detector segment (1255- 
1266 A for G130M segment B, 1330-1340 A for G130M 
side A, 1520-1533 A for G160M segment B and 1664- 
1676 A for G160M segment A). 

Both of these methods suffer from problems, especially 
for spectra with S/N ratios below ~10. In particular, due 
to the random noise, the offsets found from cross corre¬ 
lations have uncertainties of a few krn/s, so shifting lines 
to match their centroids may not yield the correct an¬ 
swer and gives rise to derived shifts that bounce around 
over short wavelength ranges. When applying the Dan¬ 


forth et al. (2014) method, in most cases the shifts are 
similar to the shifts we find using our method (described 
below), but we have found a few instances where the 
selected wavelength ranges contain a complicated set of 
IGM lines and this leads to wildly incorrect shifts - in 
one case a shift of 100 kms -1 was derived. Further¬ 
more, even after aligning in this fashion, the ISM lines 
in the combined spectra will not always be aligned with 
the Galactic 21-cm emission, as there may be general 
offsets in the spectrum chosen to provide the reference 
wavelengths. Thus, the absolute velocity scale remains 
suspect. 

To solve these problems with wavelength offsets, we 
instead take the following approach: 

(a) First, display all exposures, separately for each grat¬ 
ing and each detector side. Then, from the set of ISM 
lines usually visible in the spectrum (Fe IIAA1144, 1608 
P IIll 152, Si IIAA1190, 1193, 1260, 1304, 1526, N IAA- 
1199/1199/1200, Si III11206, S IIAA1250, 1253, 1259, OI 
11302, Ni IIAA1317, 1370,1454, C II11334, C II* 11335, 
SilV AA1393, 1402, C IVAA1548,1550, Al II11670) select 
the ones that are strong and not contaminated by IGM 
lines. We note that the N I, O I and Sill-1304 lines usu¬ 
ally are not selected because of the presence of strong 
geocoronal N I, O I and O I* emission. To this set we 
add all strong IGM lines. 

(b) Cross-correlate each line in each exposure pair, using 
a 3 A wide region. One spectrum (the first in the list) is 
chosen as a wavelength reference. 

(c) For each exposure, plot the offsets relative to the 
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flux flux 


Fig. 6.— Ratio of fitted error to CALCOS error in individ¬ 
ual 1-orbit COS exposures as function of target flux, in units of 
10 -14 erg cm -2 s —1 A - 1 , near nine ISM absorption lines. The fit¬ 
ted error is found as the rms around the polynomial fit to the 
continuum, while the CALCOS error is found by averaging the er¬ 
ror array in the original datasets over a 300 kin s' 1 wide line-free 
region near each ISM line. We include only spectra for which a 1st 
or 2nd order polynomial sufficies, and for which the continuum is 
relatively flat (variation <20% in a window several A wide). 


reference exposure as function of wavelength and fit a 
1st (usually), 2nd (sometimes) or 3rd (if the S/N ratio 
permits it) order polynomial through these points, sep¬ 
arately for side A and B of each grating. This step is 
illustrated in Fig. 5. 

(d) Apply the shift as function of wavelength to each 
spectrum and then combine (see below for more details 
on this step). 

(e) Fit the centroids of interstellar lines that are not sat¬ 
urated, not blended, not too weak and not otherwise dis¬ 
torted. Well-separated HVCs and IVCs are fitted sepa¬ 
rately. 

(f) Fit the centroid of a 21-cm emission spectrum in 
the direction of the target, with velocities referenced to 
the LSR. Comparing the centroids of the 21-cm emis¬ 
sion components with those of the ISM lines (see left¬ 
most panel in Fig. 5) then gives the wavelength offsets 
needed to align the reference spectrum. The 21-cm data 
are from the LAB survey (Kalberla et al. 2005) or the 
GASS dataset (Kalberla et al. 2010). If higher angular 
resolution data is available, these are used instead. The 
centroid determination is done separately for weak and 
strong ISM lines and separately for low- and high-velocity 
(if any) gas. For strong lines multiple components can 
blend. Usually the strong lines are not used in the end, 
since weak emission components can give substantial ab¬ 
sorption components and thus the centroids of emission 
and absorption no longer align. 

This procedure yields a final spectrum with wave¬ 
lengths such that velocities will be on the LSR scale, 


which can be converted to a heliocentric wavelength 
scale. This method works in all cases and does not re¬ 
quire the presence of STIS data. Figure 5 shows exam¬ 
ples of the offsets as function of wavelength for one of the 
targets used in this paper. 

4.2. Combining individual exposures 

To combine the aligned spectra, we add the total 
counts in each pixel and then convert back to flux, using 
the average flux/count ratio at each wavelength that was 
used in the original retrieved datafiles. This is the same 
procedure that was followed by Meiring et al. (2011) and 
Tumlinson et al. (2011), but stands in contrast to the 
procedure of Danforth et al. (2014), who instead cal¬ 
culate an inverse-variance weighted average flux. For 
cases with equal exposure time for different exposures 
and non-varying target flux, the different methods give 
the same answer. However, it is easily shown analytically 
that if the target flux varies between exposures the latter 
method gives the wrong answer. For instance, combin¬ 
ing an exposure with one in which the source brightened 
by a factor two, the inverse-variance-weights give a final 
flux in the continuum that is 4/3 times the lower value. 
But at low count rates (e.g. in the darker part of spectral 
lines), the inverse-variance-weight gives a flux more like 
3/2 times the lower value. I.e., this weighting scheme 
changes the shape of the spectral lines, which is not the 
case when combining counts. 

4.3. Error calculation 

A study of the error array in the 1-dimensional spectra 
that CALCOS provides revealed that at low fluxes the 
pipeline gives incorrect numbers. This is found by com¬ 
paring the CALCOS errors to measured errors. CAL¬ 
COS errors were calculated as the average of the error 
arrays in a 300 km s _1 wide line-free region near nine ISM 
absorption lines (in a spectrum rebinned by seven pixels). 
Measured errors are given by the rms found when fitting 
a polynomial through the line-free regions near the same 
nine ISM absorption lines, also in spectra binned to seven 
pixels. Figure 6 shows the ratio of the error array in the 
datafile to the fitted error. Clearly, for targets with fluxes 
below ~10 -14 erg cm -2 s _1 A -1 the measured rms typi¬ 
cally is lower than the value implied by CALCOS, by 
up to a factor 2 near a flux of 10“ 15 erg cm -2 s _1 A -1 . 
We have been unable to track down the source of this 
problem, but it seems possible that at low count rates 
the CALCOS pipeline overestimates one of the contribu¬ 
tions to the error calculation. 

We find a better match between calculated and mea¬ 
sured error if we just estimate the error from the Poisson 
noise implied by the total count rate. In that case the 
calculated and measured error remain similar even at the 
lowest fluxes. Our final errors are thus calculated using 
this approach, and ignore the expected, but apparently 
not actually seen, contributions from dark count, back¬ 
ground subtraction and other such items. 

5. SIMULATIONS 

5.1. Description 

In order to interpret the observations, we analyzed 
the cosmological hydrodynamical simulations of Oppen- 
heirner et al. (2010). In this section we first describe 
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R.A. 


Fig. 7. — View of the galaxies as seen from one viewpoints in the vzw simulation. Galaxies are shown by colored circles, with the same 
scalings as were used for Figs. 2 and 4. The filament axis is shown as a thick black line, found by treating the simulated galaxies in the 
same manner as the observations. Finally, an rectangular outline box shows the area of simulated sky used to measure 7V(H I) as function 
of filament impact parameter. 



p(fil) [kpc] 


Fig. 8.— Density of galaxies brighter than absolute K s =—16.3 as 
a function of filament impact parameter for each of the four simu¬ 
lated viewpoints (colors) and the observations (black). All galaxies 
out to 5 Mpc from the filament axis are included and binned every 
0.5 Mpc. The vertical scale gives the number of galaxies per square 
Mpc in each bin. 

these simulations and the information that we extracted 
from them, so that we can refer to this in the results 
section. The simulations are Gadget-2 smoothed par¬ 
ticle hydrodynamic simulations run with 384 3 gas and 
dark matter particles in a random periodic volume of 
48/i -1 Mpc. The cosmology used is Dm=0.28, Da=0.78, 
H o =70 kms -1 Mpc -1 , S2b=0.046, n=0.96, and er 8 =0.82. 


The gravitational softening length is 2.5 ft -1 kpc and the 
gas mass resolution is 3.5x 10' Mg, both of which are suf¬ 
ficient to resolve the structure of the Lyct forest. Dave et 
al. (2010) used a series of these simuluations with differ¬ 
ent feedback models to model the Lya forest statistics, 
and checked that the Lya forest statistics are resolution 
converged by using a 96 /i -1 Mpc box with 3.375x lower 
resolution. We use the “vzw” simulation, because this 
simulation’s galaxy mass function shown in Oppenheimer 
et al. (2010) is most similar to the observed galaxy mass 
function of Bell et al. (2003), which used 2MASS K s -band 
magnitudes along with SDSS magnitudes. This simula¬ 
tion mimics momentum-driven winds from starbursts as 
described in Oppenheimer & Dave (2008). However, the 
other models described in that paper give very similar 
results. 

We use SKID (Spline Kernel Interpolative Den- 
max) (http: / / www-hpcc. astro .Washington. edu/tools / 
skid.html) to identify galaxies (Keres et al. 2005; Op¬ 
penheimer et al. 2010). We determine the broadband 
photometric properties of each galaxy by summing 
up the single stellar population models of Bruzual & 
Chariot (2003) within each star particle, for which we 
have an age and metallicity, assuming a Chabrier (2003) 
initial mass function, and adding a dust correction based 
on the galaxy’s metallicity as explained in Finlator et 
al. (2006). 

5.2. Finding filaments in the simulations 

To compare these simulations to the observations we 
take a diffferent approach than usual. Usually, Simula- 
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12.8 13 13.2 13.4 13.6 13.8 14 14.2 14.4 14.6 14.8 

log N(HI) [cm - ] 


Fig. 9a. — View of the vzw simulation (with 100 kpc voxels) from four different viewpoints in the cube. H I column densities calculated 
using the Haardt & Madau (2001) extragalactic background radiation field are shown by the grey-yellow-red colors, with a contour at 
logV(H l)=13, which is our detection limit. Galaxies are shown by colored circles, with the same scalings as were used for Figs. 2, 4 and 
7. The filament axis is shown as a thick black line, found by treating the simulated galaxies in the same manner as the observations. A 
rectangular outline box shows the area of simulated sky used to measure N (H I) as function of filament impact parameter. 
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Fig. 9b.— Continued. 
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Fig. 9c.— Continued. 



Fig. 9d.— Continued. 
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Fig. 10a. — View of the vzw simulation from one viewpoint, with four different versions of the extragalactic background radiation (EGB): 
Haardt & Madau (2012) (panel a), Haardt Sz Madau (2012) times 2 (panel b), Haardt Sz Madau (2001) (panel c) and Haardt Sz Madau 
(2001) times 2 (panel d). Otherwise the meaning of the colors, contours and labels is the same as for Fig. 9. 
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Fig. 10b.— Continued. 
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Fig. 10c.— Continued. 



12.8 13 13.2 13.4 13.6 13.8 14 14.2 14.4 14.6 14.8 

log N(HI) [cm - ] 


Fig. 10d.— Continued. 
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tions are displayed by taking a 3-D cube and collapsing a 
slice onto a side of the cube. Here, however, we take the 
simulations and project the galaxies and H I volume den¬ 
sity onto the sky, taking into account the fact that a sky 
pixel represents a cone through the simulation cube, not 
a rectangular bar. We identify simulated filaments that 
look like the observed filament based on how the galaxies 
appear in the simulated skies from various locations cho¬ 
sen within the simulation, which we call “viewpoints”. 

To find viewpoints, we started with the SKID- 
identified galaxies, projecting their locations onto the 
same region of the sky as our observed filament 
(~30°x30° centered on R.A.=240°, Dec=55°). We then 
calculated each galaxy’s apparent velocity, combining its 
Hubble flow velocity with its peculiar velocity and se¬ 
lected the galaxies in a 1400 kms” 1 wide window be¬ 
tween 2900 and 4300 kms” 1 . Stepping through the sim¬ 
ulation cube in 2Mpc intervals in x, y , and 2 (i.e., 32 
steps each) and plotting the resulting distribution of sim¬ 
ulated galaxies on the sky, we found the viewpoints where 
a galaxy filament stood out visually. Next, we refined 
the search around each of these viewpoints to get the 
best viewpoint to within 1 Mpc. Using the galaxy dis¬ 
tances and simulated K s magnitudes, we find apparent 
magnitudes. Applying a magnitude cut of K s <19 and 
using just the galaxy locations, we derive filament axes 
for each viewpoint, using the same method as was used 
for the observations. 

Figure 7 presents the galaxy distribution for one view¬ 
point, in order to show that the visual impression given 
by the simulation is similar to that given by the data. 
To determine how well simulated galaxy filaments repre¬ 
sent our observed galaxy filament, we compare the galaxy 
density along the filament and the K s -band luminosity 
function (Fig. 8), using the boxes around the filaments 
seen in Fig. 9. These have lengths of 18-22 Mpc. Within 
these boxes we find 336, 157, 156 and 127 galaxies with 
K s <19 for filaments a, b, c and d, respectively. The 
lower branch of our observed filament has 137 galaxies 
and a length of 20.9 Mpc, i.e. similar to the simulated 
filaments. 

When comparing simulated and observed filament 
galaxy statistics, we apply a luminosity cut of K s =—16.3 
which corresponds to A£>=10 9 Mq and which is the 
quoted simulation resolution limit. Even though the vzw 
simulation does not quench galaxies as required by obser¬ 
vations, this does not matter, because Oppenheimer et 
al. (2010) show that the number density of galaxies above 
10 9 ° Mq is in agreement with Bell et al. (2003) within a 
factor of 30% (see their Figure 6). While we did explore 
other wind models including a no wind model and a con¬ 
stant wind model, there is statistically no difference in 
filament statistics for a given ionization background. 

Figure 8 shows the galaxy filament density as a func¬ 
tion of impact parameter for the luminosity cut. Even 
though more galaxies were used to identify the observed 
filament and make the initial identification of simulated 
filaments, we plot the galaxy densities with absolute 
K s <—16.3 binned into 500 kpc bins perpendicular to the 
filament in areal density units. The observed filament 
has a total density of 2.6 galaxies Mpc” 1 within 5 Mpc, 
while simulated filaments a), b), c), and d) have 3.8, 
2.0, 1.9, and 1.7, thus bracketing the observed filament 


density. Based on binned areal density, no filament ex¬ 
actly matches the observed filament, but filament a) is 
generally denser, and filaments b), c), and d) are less 
dense within 1.5 Mpc and have similar areal densities 
as the observations out to 4 Mpc. We also checked the 
galaxy luminosity functions along the filaments and find 
that the simulated luminosity functions are shaped like 
we expect from Oppenheimer et al. (2010). They have a 
similar density of galaxies, but more bright galaxies than 
the observations. 

5.3. Ionizing radiation in the simulations 

To calculate an H I column density map as seen from 
the selected viewpoints, we used four different levels for 
the intensity of the extragalactic background radiation. 
First, we took the model of Haardt & Madau (2001; 
HM01), in which the contribution of galaxies and quasars 
to the ionizing flux is comparable at all redshifts. Sec¬ 
ond, we used the purportedly improved model of Haardt 
& Madau (2012; HM12), which used more observational 
constraints. In this version the escape fraction of ionizing 
radiation from galaxies evolves over time such that the 
contribution of galaxies at 0 is minimal. As Kollmeier 
et al. (2014) show, the predicted H I column density dis¬ 
tribution at z=0 then fails to match the distribution ob¬ 
served by Danforth et al. (2014), giving five times too 
many H I clouds at any given column density. The 
HM01 version of the radiation field matches the observa¬ 
tions better. Kollmeier et al. (2014) discuss the possible 
causes and remedies to improve the match between the 
observed z~0 Lya column density distribution and the 
distribution implied by the HM12 model, but are un¬ 
able to come to a conclusion. Shull et al. (2015) propose 
to solve the discrepancy by using a different simulation, 
which has stronger heating and then a higher escape frac¬ 
tion (5%) of ionizing photons from starforming galaxies 
gives a factor two more ionizing photons. Khaire & Sri- 
anand (2015), however, redo the modeling of the EGB 
from scratch and suggest that the QSO luminosity func¬ 
tion alone implies a factor two more ionizing photons 
than found by HM12, which combined with a 4% escape 
fraction resolves the discrepancy. In Sect. 6.3 and 6.5 we 
describe our new way of constraining the EGB models, 
now based on the spatial distribution of the Lya forest 
in relation to galaxy filaments, rather than using the col¬ 
umn density distribution. 

We now describe our method to calculate H I col¬ 
umn densities, which depend on the assumed intensity 
of ionizing radiation. Our method produces a more re¬ 
alistic description of the observations than the standard 
method of collapsing a simulation cube onto one of the 
side planes. We start by calculating the distance to each 
(3D) simulation voxel (i.e a 100x100 kpc grid cell), given 
a viewpoint. This is converted to its Hubble flow veloc¬ 
ity, using a Hubble constant of 71 kms” 1 Mpc” 1 . To 
this recession velocity we add the peculiar velocity of the 
material in that simulation voxel to get its apparent ve¬ 
locity. Selecting a range of velocities, we then draw a 
diverging bundle of sightlines through the cube for each 
individual sky pixel, with the sightlines in a bundle close 
enough that several span the most distant voxel used. 
Sampling each sightline in steps of about 1/3 the size of 
a voxel, we sum the H I volume density that was cal¬ 
culated using value of the total gas density in that voxel 



19 


combined with the model for the EGB. This is multiplied 
by the pathlength to find the column density in each sky 
pixel. Although we use SPH simulations, we transform 
the particles onto a 3D grid with 100 kpc cells, which 
corresponds to 615 at redshifts of 3500 kms -1 . We also 
tried using 50 kpc cells, and find no statistical difference, 
which is to be expected since structures in the Lya forest 
structures are larger than this (Dave et al. 2010). 

The resulting H I column densities are shown in Fig. 9. 
From these figures it is clear that with a velocity range 
of ~1400 kms -1 a single filament as seen in the distribu¬ 
tion of galaxies is obvious, but in the H I column density 
distribution other filaments “spill over” into the selected 
velocity range. Figures 10a to d show the effect of vary¬ 
ing the intensity of the ionizing background on the H I 
column density maps. 

5.4. Properties of the gas in the simulation 

To conclude this discussion of the simulations we 
present some of the characteristics of the neutral hydro¬ 
gen. In Fig. 11 we show the relation between the total 
and neutral amounts of hydrogen for the pixels inside 
the selection box in the first filament (see Fig. 9a), using 
the HM01 EGB model. The left panel shows that within 
the about 8 Mpc wide strip along the filament the total 
hydrogen column density (integrated from cz =2900 to 
4300 kms -1 ) only varies from about log_/V(H)=18.5 to 
20.0. The concentration of dots on the right is caused by 
sightlines with high 7V(H I), which occur in voxels that 
contain a galaxy. Panel (b) reveals that the ratio of total 
to neutral hydrogen (i.e. the ionization correction) ranges 
from about 10 +65 for absorbers with logiV(H I)~13 to 
about 10 +1 when log V(H I)~19.5, with a spread of only 
about 1 clex at any given value of 7V(H I). This plot 
clearly shows that (a) the H I column density can serve 
as a proxy for the total hydrogen column density and (b) 
the H I traces only a very small fraction of the baryons, 
which was demonstrated by Dave et al. (1999), although 
they did not explicitly show this in a figure. 

The top and middle histograms in Figure 12 show the 
number of pixels and the total or neutral hydrogen mass 
as function of the total and neutral hydrogen volume 
density, for the pixels inside the outlined box for the first 
viewpoint (see Fig. 9a). The vertical line in the left pan¬ 
els represents the average hydrogen density of the uni¬ 
verse: ra(H) = 0.75x0.046 (3 H 2 )/(8 ttG) = 10 -6 ' 7 cm -3 . 
The bottom histograms give the cumulative mass frac¬ 
tion as function of volume density. Three different curves 
are shown, with the black curve including all voxels in¬ 
cluded in the outlined and the blue curve showing the 
voxels in which a simulated galaxy is present. The red 
curve is for the 27 voxels containing each galaxy (not 
double counting cases where these overlap when there are 
multiple galaxies in the list that are close together), i.e. 
they represent the circumgalactic medium out to about 
150 kpc. 

This figure shows that in a given random direction 
through the filament most voxels have low volume den¬ 
sity. The circumgalactic and galactic voxels have high 
density, as can be expected. The middle left panel shows 
that the baryon mass is fairly evenly distributed be¬ 
tween all volume densities. In the 0.1 dex bins near 
logn(H)=—8, —7, —6, —5, —4, and —3 there are 1, 3, 


6, 5, 4 and lxlO 12 solar masses of baryons, respectively. 
Thus, the total mass of baryons is dominated by the ex¬ 
tended intergalactic medium. In contrast, the distribu¬ 
tion of neutral hydrogen mass is dominated by the voxels 
around galaxies (see right middle panel). 

6. RESULTS 
6.1. Spectral analysis 

We now describe all the results we derive from the data 
and the comparison with simulations. 

We applied our wavelength-correction and error- 
calculation code to all the spectra listed in Table 1 and 
then identified the absorption lines in each, making sure 
to identify each feature between 1216 and 1240 A, i.e., 
0 to 6000 kms -1 relative to Lya. Figure 13 shows the 
relevant section of each spectrum, including a continuum 
that was determined by fitted a low-order (1st or 2nd) 
polynomial to line-free regions near the window shown. 
The targets are sorted by filament impact parameter 
(given in the bottom right corner of each panel). It is 
obvious that absorption is present close to the filament 
(left column) and absent far from the filament (right col¬ 
umn). 

Since most of our targets have relatively low redshift 
(2<0.4), there are few hard-to-identify lines and most 
features seen are clearly Lya or metal lines associated 
with a higher redshift system. In fact, there are only 
a handful of non-Lya lines that fall in the window of 
interest for this paper: Lye at 0=0.316 toward 3C 351.0, 
Ly<5 at 0=0.297 toward H 1821+643, Ly/z-Ly p in a Lyman 
limit system at ,2=0.347 toward RXJ1500.5+5517, and 
CII AA903.6, 903.9 in a Lyman limit system at 0=0.367 
toward SBS 1551+572. These lines are labeled in Fig. 13. 
All other visible features are identified as Lya. 

We used two different methods to measure the column 
densities and linewidths in the Lya lines, using the same 
continuum fits. First, we calculated column density us¬ 
ing the apparent optical depth method (see Savage & 
Sembach 1991): 

W = |iV a („)* = |^/Aln2M, 

with F(v) the observed profile, C{v) the fitted contin¬ 
uum, / the oscillator strength (0.4164 for Lya) and A 
the rest wavelength (1215.67 A). We also calculated the 
second moment of the apparent optical depth profile to 
derive a linewidth. We then deconvolved this linewidth 
for instrumental broadening (assuming an instrumental 
resolution of 20 kms -1 FWHM). Second, we used the 
VPFIT package (see Carswell et al. 2002, Kim et al. 
2007) to make a Voigt profile fit to the Lya lines. VPFIT 
version 10.2 (http:// www.ast.cam.ac.uk/ffc/vpfit.html) 
was used. The theoretically calculated line spread func¬ 
tion for the COS G130M grating at Lifetime Position 
1 (Kriss 2011) was used to convolve with the model fit 
profile. This line-spread function accounts for both scat¬ 
tering in the far wings due to microroughness in the pri¬ 
mary mirror and zonal polishing errors in the primary 
and secondary mirrors, which results in a COS LSF with 
reduced core intensity and non-Gaussian strong wings. 

Equivalent width detection limits for non-detections 
were found as three times the error in the equivalent 
width for a line that is 50 kms -1 wide. This was then 
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Fig. 11.— Scatter plot comparing the total baryon column density (log 1V(H)) against the column density of neutral hydrogen (log .V( T1 I)) 
for the rectangular box outlining the filament in Fig. 9a. I.e., the column density when integrating from cz =2900 to 4300 kms" 1 (a 20 Mpc 
pathlength), which is dominated by the voxels inside the filament. The concentration of points around log IV (H l)~18 in the left panel 
corresponds to pixels with a galaxy present. The panel on the right shows how the fraction of hydrogen that is neutral varies with the 
observed H I column density; i.e. it represents the inverse of the ionization correction. 
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Fig. 12.— Top panels: number of 100 kpc 3 voxels as function of total or neutral hydrogen volume density in the rectangular box outlining 
the filament in Fig. 9a. Middle panels: total mass of hydrogen in each 0.1 dex wide bin in logn(H) (left) or logn(H) (right). Bottom 
panels: cumulative mass fraction as function of volume density. The black curves are for the case that includes all voxels in the box. The 
blue curves are for voxels in which SKID identified a galaxy, while the red curves are for the 27 voxels around those galaxy-containing 
voxels (i.e. a 300x300 kpc box). The vertical line in the left two panels indicates the average density of hydrogen (logn(H)=—6.6, while 
the vertical line in the right panels indicates the density that corresponds to a H I column density of 10 13 cm -2 in the 20 Mpc deep box. 


converted to a column density detection limit by calculat¬ 
ing the N a (v) integral of a line with the equivalent width 
equal to the detection limit and an FWHM of 50 kms -1 . 
On average this limit is about 10 13 cm -2 . 

Table 3 presents the measured equivalent widths, col¬ 
umn densities and linewidths determined using both 
the apparent optical depth and profile fitting meth¬ 
ods. For the sake of completeness, the table includes 
Lya detections for the extended velocity range shown 


in Fig. 13 (2400 kms -1 to 4800 kms -1 ). In this figure 
the Lya components that we associate with the galaxy 
filament are shown by blue horizontal lines extending 
from v(Lya)-2b(Lya) to u(Lya)+26(Lya). Four com¬ 
ponents can be considered Broad Lya Absorbers (BLAs; 
6>40 kms -1 ; see Richter et al. 2004; Lehner et al. 2007). 
The four components in Table 3 that are not in the ve¬ 
locity range where filament galaxies are found are shown 
by green horizontal bars. 
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Fig. 13. — Relevant sections of the COS spectrum for each of the 24 targets analyzed. The targets are sorted by filament impact parameter, 
given in units of kpc in the bottom right corner. All absorption features are Lya, except for the five labeled Ly5, Lye, Ly£, Ly-lim and C III. 
The blue dashed vertical lines shows a 1200 kms —1 wide window around the velocity of the nearest filament axis segment (see Sect. 2). 
Only lines falling inside this window are considered to be associated with the filament. These are indicated with blue horizontal bars above 
the Lya absorption. Lya lines that are not considered associated with the filament are indicated by green horizontal bars. The four BLAs 
are indicated by the thicker blue bars. 
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TABLE 3 

Ly a PARAMETERS AND DETECTION LIMITS FOR ABSORPTION BETWEEN 2400 AND 4800 KM S —1 


Target 

(1) 

Flux 

[f.u.] 

(2) 

S/N 

(3) 

p(fil) 

kpc 

(4) 

v(Lya) 

(N a (v)) 

[km/s] 

(5) 

EW 

[mA] 

(6) 

log 7V(H I) 
(JV„(w)) 

[cm 2 ] 

(7) 

(Na(v)) 

[kms x ] 
(8) 

v(Lya) 

fit 

kms -1 

(9) 

log AT(H I) 
(fit) 
[cm -2 ] 

(10) 

b 

(fit) 

[kms -1 ] 

(11) 

Note 

(12) 

3C309.1 

0.12 

13.3 

3908 


<39 

<12.85 






3C351.0 

1.28 

8.4 

569 

3596+2 

185+12 

13.68+0.03 

33.1+2.2 

3597+3 

13.77+0.05 

28.9+4.0 

b 





3456+3 

140+15 

13.49+0.05 

44.6+2.5 

3459+5 

13.53+0.06 

39.9+8.0 


4C63.22 

0.33 

8.2 

2333 

2421+2 

228+26 

13.83+0.07 

35.7+2.7 

2420+4 

13.90+0.08 

32.6+6.3 

a 

FBS1526+659 

0.59 

12.0 

369 

3462+4 

136+19 

13.48+0.06 

46.4+2.3 

3476+7 

13.48+0.08 

38.6+9.7 


H1821+643 

3.72 

52.4 

3249 

2824+4 

44+4 

12.93+0.04 

65.0+2.2 

2825+6 

12.96+0.04 

62.5+8.4 

a 





4087+2 

36+3 

12.86+0.03 

19.4+2.2 

4087+1 

12.91+0.03 

18.7+2.6 


Kaz447 

0.30 

11.5 

184 


<48 

<12.94 






Mrk290 

2.58 

30.4 

97 

3085+2 

505+6 

14.29+0.01 

56.9+2.9 

3089+1 

14.37+0.01 

52.9+1.5 






3202+2 

318+5 

14.04+0.01 

37.8+3.4 

3204+1 

14.07+0.02 

32.1+1.4 


Mrk486 

0.18 

10.0 

1910 

4387+3 

161+17 

13.63+0.06 

25.4+1.4 

4386+3 

13.74+0.06 

25.0+4.7 

a 

Mrk817 

10.83 

47.1 

6208 


<9 

<12.19 






Mrk876 

4.00 

48.8 

540 

3472+2 

280+4 

13.90+0.01 

53.2+3.9 

3476+0 

13.92+0.02 

24.7+0.9 

c 









3470+5 

13.22+0.07 

80.9+11.1 


PG1626+554 

2.66 

28.6 

5684 


<21 

<12.52 






RBS1483 

0.76 

11.2 

2260 

2721+2 

259+19 

13.85+0.04 

44.6+3.1 

2726+2 

13.90+0.03 

35.0+2.9 

a,d 

RBS1503 

0.78 

13.2 

306 

3272+2 

270+16 

13.88+0.03 

44.1+2.2 

3269+3 

13.81+0.08 

26.9+5.3 

e 





3388+2 

56+15 

13.06+0.10 

39.9+3.0 

3306+23 

13.59+0.12 

116.8+28.8 


RXJ1500.5+5517 

0.22 

13.8 

2031 

3595+2 

120+13 

13.44+0.05 

28.0+2.7 

3592+4 

13.51+0.06 

26.8+6.3 

f 

RX J1503.2+6810 

0.75 

12.7 

1016 


<51 

<12.96 






RXJ1508.8+6814 

0.24 

8.0 

1202 


<102 

<13.27 






RXJ1608.3+6018 

0.81 

20.2 

880 

2983+4 

388+8 

14.15+0.02 

43.0+2.6 

2983+1 

14.26+0.03 

37.0+2.2 






2877+2 

117+10 

13.42+0.04 

37.7+6.0 

2886+5 

13.46+0.06 

37.4+7.0 


RXJ1717.5+6559 

0.23 

9.6 

4330 

4705+2 

234+16 

13.86+0.04 

29.6+3.3 

4705+2 

13.86+0.03 

32.3+2.4 

a 

SBS1458+535 

0.25 

11.4 

563 


<57 

<12.99 






SBS1503+570 

0.31 

13.5 

2395 


<51 

<12.94 






SBS1521+598 

0.18 

10.0 

2390 


<66 

<13.08 






SBS1537+577 

0.14 

8.8 

293 

3541+2 

436+26 

14.20+0.05 

49.4+3.0 

3541+3 

14.43+0.18 

37.6+6.1 






3260+2 

360+34 

13.98+0.05 

70.9+3.9 

3257+7 

14.03+0.05 

73.2+9.8 


SBS1551+572 

0.11 

7.3 

1778 

4097+3 

195+21 

13.74+0.07 

28.3+1.6 

4097+5 

13.93+0.10 

30.9+7.6 


SBS1624+575 

0.31 

12.4 

4103 


<48 

<12.94 







Note. — Col. 1: Target name; Col. 2: Flux in units of 10 -14 erg cm -2 s -1 mA -1 ; Col. 3: Signal to noise ratio at C2=3000 kms -1 ; Col. 4: cz 
of Lya detection as found from N a (v ) integral; Col. 5: Equivalent width of Lya detection; Col. 6; Column density of Lya detection found as the 
integral of the N a (v ) profile; Col. 7: b -value of Lyo: detection found from the second moment of the N a (v ) profile; Col. 8: cz of Lya detection as 
found from profile fitting; Col. 9: Column density of Lya detection as found from profile fitting; Col. 10: 6-value of Lya detection as found from 
profile fitting; Notes: (a) Component not considered to be associated with the filament; (b) The 3C351.0 spectrum was taken using the STIS-E140M 
echelle grating; (c) The very high S/N spectrum of Mrk876 is best fitted using both a narrow and a broad component, but with the AOD method 
only a single measurement can be made; (d) This line may have a second component at lower velocities, although a two component fit is not easily 
justified - yet the AOD 2nd moment is increased; (e) The properties of the second component are uncertain because it is weak and blended with 
the stronger component, see discussion in the text; (f) This line lies 0.6 A above the Lyman limit of an absorber with logiV(H i)~16.5, i.e. at a 
wavelength where the individual Lyman lines are blended and weak. 
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Fig. 14.— Comparison of linewidth (left) and column density (right) derived using the apparent optical depth (AOD) method (horizontal 
axes) and Voigt profile fitting (vertical axes). Closed symbols are for single-component Lya absorption lines, while open circles are for 
individual components in multi-component lines. 


Toward Mrk290 (filament impact parameter 97 kpc) 
two components are clearly visible. Associated O VI is 
also seen in the FUSE spectrum of this target. These 
components were analyzed in detail by Narayanan et al. 
(2010), who concluded that they originate from a cloud 
with iV(H)=4xl0 19 cm -2 at T=1.4xl0 5 K which is also 
photoionized. Where Narayanan et al. (2010) focused on 
a possible association with NGC 5987 (impact parameter 
475 kpc (1.9 R v ir see Table 2), it is quite likely that in 
fact this absorber provides evidence for the presence of 
WHIM gas in the filament. 

The spectrum of Kaz447 (filament impact parameter 
164 kpc) appears to show a small dip near 4300 kms -1 . 
This feature is not significant, however, measuring as 2er. 

For SBS 1537+577 our S/N ratio is insufficient to be 
sure that the apparent broad component can only be fit¬ 
ted with a single gaussian or whether it is a combination 
of two narrower gaussians. We proceed assuming it is a 
single broad component. 

Toward RBS1503 two components are clearly 
present. The profile fitting method gives 6-values of 
26.9+5.3 kms -1 and 116.8+28.8 kms -1 , but the two 
components are heavily blended and the S/N ratio of the 
spectrum is not sufficient to make a reliable fit. There¬ 
fore, in Table 3 we list a very different velocity for the 
second component in Col. 4 vs Col. 8. This is due to 
the fact that when using the N a (v) method we can only 
measure the line wing, which has a centroid velocity of 
3388 kms -1 . 

The BLA component toward Mrk876 overlaps with a 
narrower component at the same velocity, but because 
the S/N ratio of the spectrum is very high (~50), it is 
clearly significant and not due to the wings on the COS 
line spread function. 

It is clear from Fig. 13 that the strongest and widest 
absorbers occur toward sightlines passing closest to the 
filament axis (within 660 kpc), as is the case for multi- 
component absorbers. The exceptions are a weak line to¬ 
ward FBS 1526+659 and non-detections toward Kaz 447 
and SBS 1458+535. The former maybe due to the fact 


that FBS 1526+659 is at the far end of the upper fil¬ 
ament and the galaxy density near it is rather low. So 
even though it is nominally close to the filament axis, it is 
not in a region with a lot of gas. In the other two cases 
we must be looking through a hole, with the Kaz 447 
non-detection being especially interesting and revealing 
small-scale structure because it is close to the siglitline 
to 3C 351.0 where two components are seen. 

A map of the galaxies between 1500 and 2900 kms -1 
shows that the non-filament absorbers toward 4C 63.22, 
RBS 1483, and H1821+643 can be associated with an¬ 
other structure, which is clearly visible, but not quite as 
well defined as the filament shown in Fig. 4. Similarly, 
the component at 4705 kms -1 toward RX J1717.5+6559 
is at the end of yet another filament that extends from 
(R.A, Dec.)~(260°, 65°) to ~(240°, 45°) with cz between 
5000 and 6200 kms -1 . Only the absorber at 4387 kms -1 
toward Mrk486 is an orphan that is not clearly associ¬ 
ated with any filaments, although it is not too far off from 
our main filament (off by 1.9 Mpc and 786 kms -1 ) and 
a case could be made that it just represents an outlier 
and should be included. 

We note that all our detections have H I column den¬ 
sities below about logAl(H I) =14.5, whereas absorbers 
that are associated with galaxy halos typically have 
higher column densities. 

In Fig. 14 we compare the column densities and 
linewidths derived using both the apparent optical depth 
and profile fitting methods. This reveals that the decon¬ 
volved apparent-optical-depth linewidths are typically 
wider than those derived from profile fitting - the ra¬ 
tio 6(AOD)/6(fit) is 1.08+0.11. On the other hand, 
column densities measured using the apparent optical 
depth method are consistently lower than those found 
from profile fitting, with AlogAl=0.07±0.06 dex. Both of 
these discrepancies can be understood as a consequence 
of the non-gaussian shape of the COS linespread function 
(LSF). Since the LSF has strong line wings, the appar¬ 
ent line width becomes larger than what it would be if 
the LSF was a single gaussian, even after deconvolution. 
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Similarly, the LSF smearing increases the apparent flux 
in the center of the line, lowering the apparent optical 
depth. As our comparison shows, this ends up as an 
about 10% increase in the apparent linewidth and an 
about 0.1 dex decrease in the derived column density. 

6.2. Equivalent width and line width vs filament impact 

parameter 

Using the galaxy and filament impact parameters given 
in Table 2, we can now plot the distribution of Lya lines 
relative to galaxies and relative to the filament axis. This 
is done in Fig. 15. We use the equivalent width instead 
of the column density of the Lya absorption to facilitate 
a comparison of the results for the 24 sightlines in this 
paper with the much larger sample of 125 Lya detections 
given in Wakker & Savage (2009). 

The two top panels (a,b) of Fig. 15 correspond to the 
conventional way of analyzing the relation between Lya 
and galaxies. For non-detections we find the galaxy im¬ 
pact parameter as that to the nearest galaxy with veloc¬ 
ity between 2900 and 4300 kms” 1 , i.e. galaxies shown 
in Fig. 4. For Lya absorbers the relevant impact pa¬ 
rameter is that to the nearest galaxy whose velocity is 
within ±400 kms” 1 of that of the absorber. As Ta¬ 
ble 2 shows, there are only a few cases where the nearest 
galaxy is very small (L<0.1 L*; D< 7.2 kpc) (Mrk486, 
SBS 1503+570, SBS 1521±598, SBS1624+575), but in 
all of those cases the impact parameter to the nearest 
galaxy with L>0.1L* is not very different, so we avoid 
the possibly confounding problem of whether to choose 
the small impact parameter to a dwarf galaxy or the 
much larger value to a substantial galaxy. 

In this figure, we use a red symbol if the sightline 
passes within 1.5 virial radii of a galaxy (which happens 
in just one case, for SBS 1537+577), and a blue symbol 
otherwise. There are four sightlines (3C 351.0, Mrk876, 
RXJ1608.3+6018 and SBS 1551+572) where the nearest 
galaxy is between 150 and 250 kpc, but the impact pa¬ 
rameter also is between 1.7 and 2.0 virial radii, so it is not 
likely that the absorption is associated with the galaxy 
halo. 

Panels (a) and (b) of Fig. 15 show that the sample 
of sightlines in this paper has similar properties as the 
larger sample in Wakker & Savage (2009). I.e., (1) the av¬ 
erage equivalent increases for smaller impact parameters; 

(2) every sightline with a galaxy impact parameter below 
250 kpc shows a Lya line; (3) for some high equivalent 
width absorbers the nearest galaxy is separated by more 
than 250 kpc; (4) non-detections become more common 
at larger impact parameters; (5) the typical linewidth 
is larger at smaller impact parameters. Detections at 
large impact parameters are typically noted, but their 
properties are mostly ignored because they can’t easily 
be interpreted in terms of galaxy halos. Such absorbers 
only play a role when interpreting the column density 
distribution of Lya lines in order to estimate the baryon 
content of the Lya forest. 

Panels (c) and (d) give a different way of looking at 
these Lya absorbers and non-detections. Here the red 
symbols indicate the components toward SBS 1537+577 
which originate at 1.0 times the virial radius from a 
galaxy and thus conceivable could be associated with 
that galaxy. The closed blue stars show the components 
in the lower filament branch seen in Fig. 4, while the 


open blue stars are for the upper branch. It is clear that 

(1) The high equivalent width absorbers occur close to 
the filament axis, whether or not there is a galaxy within 
one virial radius of the sightline. 

(2) There are no detections with equivalent width 
>50 mA further than 2 Mpc from the filament axis, un¬ 
like what is the case in panel (a) where for some absorbers 
with equivalent width >100 mA there is no galaxy within 
2 Mpc. 

(3) There is a strong correlation between equivalent 
width and filament impact parameter. This is especially 
obvious for the lower branch (closed blue stars). Only 
the broad component toward RBS 1503 (p(fil)=373 kpc, 
EW=56 mA) does not follow this pattern, but its equiv¬ 
alent width cannot be measured properly and the 
linewidth is very uncertain. 

(4) The increase in linewidth with decreasing filament 
impact parameter is more systematic than the increase 
in linewidth with decreasing galaxy impact parameter. 

(5) The four Lya absorbers that fit the definition of 
a BLA (6>40 kms” 1 ) all occur in sightlines with two 
components and they only occur at the smallest filament 
impact parameters. This indicates that there the struc¬ 
ture of the gas gets more complicated and/or that the 
gas gets hotter. With just four cases out of twelve sight¬ 
lines we will need more sightlines and more filaments to 
make this a statistically sound conclusion, however. 

6.3. N(H I) vs filament axis: observations vs 
simulations 

To compare observed H I column densities against the 
simulations, we first show a scatter plot of N(H I) against 
filament impact parameter, for four viewpoints for the 
vzw model (Fig. 16) using the HM01 version of the EGB. 
The four viewpoints are sorted by the density of galaxies 
along the filament axis (3.8, 2.0, 1.9 and 1.7 galaxies per 
Mpc for panels a, b, c and d, respectively). Compare to 
2.6 galaxies per Mpc for the observed filament. In Fig. 17 
we show the effect of varying the EGB for one viewpoint 
for the vzw model. Observed values are included as col¬ 
ored stars, using the same coding as in Fig. 15. 

The orange points in Fig. 16 represent H I column den¬ 
sities that lie above the average detection limit of the ob¬ 
servations (10 13 cm” 2 , see Sect. 6.1), for directions that 
lie inside the rectangular box outline in Fig. 9. The black 
lines give the 10th, 25tlr, 50th, 75th and 90th percentile 
of column density. The comparison between the differ¬ 
ent viewpoints shows that there are differences in detail 
between the four simulated filaments, but in general the 
column density is always above the detection limit at low 
filament impact parameters (<~500 kpc, while it drops 
at higher impact parameters. For three of the viewpoints 
almost 100% of the sightlines have 1V(H I)> 10 13 cm” 2 at 
zero impact parameter. Across the filament the 10th to 
90th percentile of column densities have a spread of about 
a factor of 100. For three viewpoints column densities 
>10 13 cm” 2 correspond to the 50% or 75% line at sev¬ 
eral Mpc from the axis, while for one viewpoint the 75% 
percentile is crossed at 2 Mpc. The observed sightlines 
appear to provide a fair sample of the expected column 
densities. 

We note here that we adopted the momemtum-driven 
wind (‘vzw”) version of the simulation of Oppenheimer 
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Fig. 15.— Ly a equivalent width (left) linewidth (right) vs filament impact parameter to galaxies (top) and to the filament (bottom). 
Five sight lines have two components, which are shown separately. For the top panels, the impact parameter is that to the nearest galaxy of 
any size with velocity within ±500 kms" 1 of a Lya detection. Open black circles show the results from Wakker & Savage (2009). The red 
star is for the one sightline that passes within one virial radius from a galaxy. The closed blue stars are for sightlines in the “lower branch”, 
running from (RA,Dec)=(220,47) to (247,64), while open blue stars are for the “upper branch”, running from (270,60) to (230,68). Blue 
downward pointing triangles are upper limits for sightlines with no detected Lyo; in the velocity range 2900 to 4300 kms -1 . The vertical 
lines in panels (a) and (b) are at an impact parameter of 250 kpc. The horizontal line in the left panels indicates the typical equivalent 
width detection limit of 50mA, while the one in the right panels show the canonical 6=40 kms -1 separation between narrow and broad 
lines. Two non-detections from Table 3 are missing in the bottom plot since they have filament impact parameter >5 Mpc. Note that the 
detection at p(fil)=3300 kpc does not count for the statistic shown in Fig. 18 as it has log7V(H l)<13. The trend of increasing equivalent 
and line width with decreasing filament impact parameter is clear. The exceptions (the points at 50 mA and 28 kms -1 at 400 kpc) come 
from the multi-component feature toward RBS 1503 that is difficult to fit. Note that for all but four of the blue points the nearest galaxy 
is at least 500 kpc distant, or more than 3R v ir- For three the nearest galaxy is at 200-350 kpc (1.5-2.5 R v ir) and the red point is at 91 kpc 

(l.ORvir). 


& Dave (2008) as a reference, rather than the no-wind, 
constant-wind or energy-driven wind version. However, 
although in these simulations the galactic winds heat and 
disturb the IGM close to galaxies, at the Mpc scales of 
the galaxy filament they have little influence. We did 
make the equivalent version of Fig. 16 using different 
wind models, but in all cases the distribution of V(H I) vs 
filament impact parameter is essentially the same, vary¬ 
ing only by at most 10% in any individual sightline. 

In Fig. 17 we analyze the impact of the assumed inten¬ 
sity of the extragalactic background radiation showing 
the scatter plot of H I column density vs filament impact 
parameter from the viewpoint shown in Fig. 9a and the 
top left in Fig. 16 (i.e. the filament with the best match 
between observed and simulated galaxy density). We use 
four different versions of the Haardt & Madau model: 
HM01 and HM12, and both scaled up by a factor two. 
As expected, the assumed ionizing radiation field has a 


huge impact on the predicted column densities, with the 
HM12 version predicting almost no column densities be¬ 
low 10 13 cm -2 and the HM01x2 version predicting thin 
filaments at the 10 13 cm -2 level. 

6.4. Filament gas vs galaxy mass 

Combining our measurements of H I column density 
with the ionization corrections implied by the simulation, 
we can estimate the total mass of baryons in the fila¬ 
ment. Fitting the relation between ionization correction 
and H I column density that is shown in the right panel of 
Fig. 11 yields log IV(H) = 0.36 log iV(H I) + 14.37, with 
a typical variation on the order of ±0.3 dex. Note that 
this relation is implied when using the HM01 version of 
the EGB. The coefficient will change in proportion to 
the strength of the ionization radiation field. Using this 
conversion, the observed H I column densities listed in 
Table 3 imply total hydrogen column densities ranging 
from logV(H)=19.2 to 19.6, with an average of 19.35. 
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Fig. 16. — Scatter diagram of H I column density vs filament impact parameter for the four different viewpoints in the vzw model shown in 
Fig. 9. The nominal Haardt & Madau (2001) EGB prescription was used. Only points inside the rectangular outline box in that figure are 
shown. Column densities above 10 13 cm -2 are shown by orange points, lower column densities by yellow points. The downward pointing 
triangles give observed upper limits as function of filament impact parameter, while stars show the detections, with the red star for the 
sightline within 150 kpc of a galaxy. The black lines give the 10th, 25th, 50th, 75th and 90th percentile of column densities in intervals of 
200 kpc. 
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Fig. 17.— Same as Fig. 16, but comparing the vzw model from one viewpoint but with five different versions of the intensity of the EGB, 
as labeled. 
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Fig. 18.— Detection fraction as function of filament impact parameter for different four viewpoints inside the simulation cube. Differently 
colored lines indicate four different versions of the EGB, as shown by the label in the bottom left panel. The detection fraction is the 
fraction of sightlines in a 200 kpc (for simulations) or 500 kpc (for observations) interval where 7V(H l)>10 13 cm -2 . The observed detection 
fractions (black line) are 80% (4/5) for p=0-0.5 Mpc, 75% (3/4) for p=0.5-l Mpc, 0% (0/2) for p= 1—1.5 Mpc, 50% (1/2) for p= 1.5—2.0 Mpc, 
and 25% (1/4) for p= 2-2.5 Mpc. No detections are found for seven sightlines with p= 2.1—5.0 Mpc, as indicated by the thick line at zero. 
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We note that for the absorber toward Mrk 290 Narayanan 
et al. (2010) used the combination of O VI and Lya and 
a hybrid ionization model to derive logIV(H)=19.14. For 
this cloud the total hydrogen column density implied by 
the simulation is log 7V(H)=19.44, which is within the 
±0.3 dex range implied by the log?r(H) vs log?r(H I) cor¬ 
relation. If the ionizing background were twice as strong, 
the two values would match. 

In Sect. 3.2 we estimated the total area of the filament 
as 170 Mpc 2 . With the average loglV(H)=19.35 this then 
implies a total hydrogen mass in the gaseous filament of 
3x 10 13 M q , or 5.2x10 13 Mq in baryons. 

We can compare this to the mass of the galaxies in 
the filament. This can be derived by using the simulated 
galaxies, from which we find that log (mass in Mq) = 
log (luminosity in L„) + 11.1. Combining this with the 
galaxy luminosities derived in Sect. 2.1 and adding the 
implied masses for all galaxies inside the stripe boxes 
shown in Fig. 4 gives a total baryonic galaxy mass of 
1.4x 10 13 Mq. 

Thus, although we use the simulations to derive the 
conversions between luminosity and mass and between 
neutral and total hydrogen, scaling the observed galaxy 
luminosities and H I column densities implies that the 
21 absorbers in our sample represent four times more 
baryons than all the baryons inside all the galaxies that 
were used to define the filament. We previously found 
that half of the Lya absorbers are within 400 kpc of a 
galaxy (see Table 9 in Wakker & Savage 2009), implying 
that the baryonic mass of the circumgalactic medium 
is comparable to that of the intergalactic gas, and that 
both are several times more massive than the condensed 
baryons inside galaxies. 

6.5. Detection fraction 

With the limited number of sightlines we have avail¬ 
able at the moment, we cannot make the full distribu¬ 
tion of column densities in different impact parameter 
intervals. But we can compare the observed and pre¬ 
dicted detection fraction as function of filament impact 
parameter. This is shown in Fig. 18. To make this fig¬ 
ure we calculated the fraction of sightlines for which the 
vzw simulation predicts an H I column density above 
10 13 cm -2 in intervals of 200 kpc, using the Haardt & 
Madau (2001, 2012) EGB model. We also calculated 
the fraction of sightlines for which we detected H I, in 
intervals of 500 kpc. The 10 13 cm -2 boundary value is 
the average detection limit for our data (see Sect. 6.1). 
In all but two sightlines we are sensitive to Lya with 
logIV(H I) above this limit. There are only six sight¬ 
lines where the S/N ratio is >20 and the detection limit 
is lower than log IV(H I)=12.9, (H1821±643, Mrk 290, 
Mrk817, Mrk876, PG 1626±554 and RX J1608.3±6018). 
In this set there is just one H I line with log IV(H I) below 
the limit (12.9 toward H 1821±643, see Table 3). For the 
statistics discussed here we discount this line as it has 
logIV(H I) below the 13.0 limit. 

The observed detection rate (black histogram in 
Fig. 18) is based on relatively few sightlines, and there¬ 
fore the detection fraction as function of filament impact 
parameter remains fairly uncertain. See the figure cap¬ 
tion for the percentages and numbers of sightlines. The 
three robust conclusions that can be drawn are that: 

(1) At filament impact parameters below 1 Mpc the de¬ 


tection rate is high: seven out of nine sightlines show 
Lya. 

(2) The detection rate regularly decreases with filament 
impact parameter. 

(3) No Lycc lines are found in the seven sightlines with 
impact parameter >2.1 Mpc. 

Comparing the observed detection fraction to the pre¬ 
dicted fraction shows that the HM12 model (blue line) 
grossly overpredicts the number of Lya absorbers, i.e. the 
IGM is underionized in this model. Even scaling it up 
by a factor two (magenta line), as proposed by Shull et 
al. (2015) does not help. The orange line (corresponding 
to HM01, which has a factor 3.7 more ionizing photons) 
corresponds much more closely to the observed detec¬ 
tion rate at small filament impact parameters (<2 Mpc). 
However, even then the detection fraction along the axis 
is predicted to higher than the observed 80%. Further, 
some Lya lines are expected at filament impact parame¬ 
ters >2 Mpc. For the best matched filaments (a) and (b), 
the expected probability given the orange and red lines 
in Fig. 18 ranges from ~20% to ~40%. Given this expec¬ 
tation, the probability to find zero detections in a sample 
of seven ranges from (1—0.4) 7 to (1—0.2) 7 or ~3%-20%. 
Thus, the fact that we do not find any Lya lines at large 
impact parameter has a low probability. 

In summary: the HM12 version of the EGB has too few 
ionizing photons, while the shape of the run of detection 
fraction with filament impact parameter may not match 
the observations. 

7. DISCUSSION 

From our results it is clear that understanding the ma¬ 
jority of Lya absorbers is best done in terms of large- 
scale structure, rather than in terms of galaxy halos. Al¬ 
though strong (IV(H I)>10 15 cn± 2 ) absorbers are asso¬ 
ciated with the circumgalactic medium of galaxies, the 
weaker Lya lines are more likely to be associated with fil¬ 
aments. Studies of the galaxy-galaxy vs galaxy-absorber 
correlation function (e.g. Tejos et al. 2014 and references 
therein) also indicate that the two differ. As can be seen 
from Table 2, only one of our sightlines passes within 
one virial radius of a galaxy (SBS 1537±577 with impact 
parameter 91 kpc to SDSS J153802.75±573018.3, which 
has R v ; r =87 kpc according to our approximation to the 
relation of Stocke et al. (2013). In general, we previ¬ 
ously found that 50% of Lya absorbers occur more than 
300 kpc away from the nearest galaxy. As Fig. 15 shows, 
however, the properties of Lya absorbers far from galax¬ 
ies are not random, but they correlate well with the im¬ 
pact parameter to the galaxy filament that we studied. 
Not only are the largest equivalent widths seen closest to 
the filament axis, there is also a hint that the linewidth 
decreases with filament impact parameter. However, our 
sample is not large enough to come to a definitive con¬ 
clusion on this. 

With our relatively small sample, the cleanest way to 
compare observational results to simulations is in terms 
of the detection fraction. As Fig. 18 shows, the fraction of 
sightlines that show Lya absorption is maximal near the 
filament axis and decreases rapidly away from the axis, 
while no detections are found more than 2.1 Mpc away 
from the filament axis. Further, the expected decrease 
in detection fraction depends strongly on the intensity 
of the ionizing flux in the extragalactic background ra- 
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diation field (EGB). Using the most sophisticated model 
that fits the high redshift Ly a forest (Haardt & Madau 
2012) yields an ionizing flux that is far too low at z= 0, 
by a factor four to five. This conclusion was previously 
reached by Kollmeier et al. (2014), who compared the ob¬ 
served column density distribution (from Danforth et al. 
2014) with the predictions. Where Shull et al. (2015) find 
a factor two more intense ionizing flux by using a differ¬ 
ent hydrodynamical simulation and suggesting a higher 
escape fraction from galaxies for Lyman continuum pho¬ 
tons, Kliaire & Srianand (2015) suggest that the contri¬ 
bution from QSOs is twice as high as Haardt & Madau 
(2012) had found and suggest a 4% Lyman continuum 
escape fraction to imply an ionizing flux at z =0 that is 
five times higher than HM12. 

Our results show that the problem is even more se¬ 
vere, in that the models predict the wrong shape for the 
run of detection fraction with filament axis. Matching 
the detection fraction at low impact parameters means 
that there is too much H I at filament impact parame¬ 
ters >2 Mpc for simulated filaments (a) and (b), which 
have galaxy density similar to that in the observed fil¬ 
ament. Matching to the non-detections at large impact 
parameters means that too few absorbers are predicted 
at low impact parameters (see e.g. filament (d) in Fig. 18. 
This could indicate that ionizing radiation field is even 
stronger than expected far from filament axis, or that 
the total volume density of hydrogen is lower than the 
simulations predict. 

To follow up on this work, we hope to collect other 
simulations, using different prescriptions for galaxy for¬ 
mation and possibly different predictions for the distri¬ 
bution of the gas around filaments. We are also working 
on increasing our observational sample by analyzing the 
nine other nearby galaxy filaments that we have located. 
None of those filaments is sampled by more than 10 sight¬ 
lines, but cumulatively we will be able to more than dou¬ 
ble our sample, using COS archival data. It would also 
be useful to observe of one these filaments with a much 
denser pattern of sightlines. 

8. SUMMARY 

We present a study of a local filament in the Lya for¬ 
est in three dimensions. We obtained spectra of 24 AGN 
obtained with the Hubble Space Telescope, and we de¬ 
scribe a new method to properly align multiple expo¬ 
sures taken using the COS instrument. We measure the 
properties of 21 Lya absorbers with cz between 2400 and 
4800 kms -1 , seen in 17 of the sightlines. We associate 15 
absorbers with a filament of galaxies that is apparent in 
the distribution on the sky for galaxies with cz =2900 to 
4300 kms -1 . We present a method to objectively derive 
an axis for the filament and analyze the Lya absorbers 
with reference to this axis. We also search for similar 
galaxy filaments in an SPH simulation and determine fil¬ 
ament axes in the same manner as was done for the real 
sky. We then compare the observed and predicted H I 
column density as function of the sightline and trans¬ 
verse dimensions, for several different prescriptions for 
the strength of the intergalactic ionizing radiation field. 
We derive the following conclusions. 

8.1. Data Handling 


(1) We describe a method to correct the COS wave¬ 
length scale that aligns individual exposures without as¬ 
suming a constant offset, by crosscorrelating all ISM and 
IGM lines in each exposure. Further, in contrast to other 
methods that have been described, we properly deter¬ 
mine an absolute wavelength scale for the combined spec¬ 
trum by aligning the ISM lines with a 21-crn spectrum. 
This prevents smearing of absorption lines in misaligned 
spectra and allows us to properly assess the alignment of 
lines in intergalactic absorption-line systems. 

(2) An analysis of the error array produced by CAL- 
COS shows that those errors tend to be higher than 
the errors measured from the rms around a fit to the 
continuum in the spectra. For target fluxes below 
~10 14 erg cm -2 s -1 A -1 this results in an overestimate of 
the errors - on average by 50% for the faintest targets. 
Thus, using CALCOS errors can lead to large discrepan¬ 
cies when determining the significance of weak lines and 
of detection limits. 

(3) Using the apparent optical depth method to mea¬ 
sure column densities and linewidths in COS spectra 
leads to values that are 10% lower and 10% higher, re¬ 
spectively, when compared to measuring these quantities 
using profile fitting; the latter more properly takes into 
account the complex COS line-spread function. 

8.2. Properties of Lya lines 

(4) Our set of sightlines does not generally sample 
galaxy halos. All but one of the Lya lines originate 
far from galaxies 13 out of 15 absorbers (in 10 sight¬ 
lines) in our sample do not pass within 1.5 virial radii 
of any galaxy with luminosity greater than our com¬ 
pleteness limit of 0.05 L* (about the luminosity of the 
SMC). We use 1.5 virial radii as the criterion since this 
is about the radius at which theoretical considerations 
(Oort 1970, Mailer & Bullock 2004) place the boundary 
between infalling and intergalactic gas. A fair number of 
fainter galaxies are also included in our sample, but none 
is close to a sightline. Using a simple impact parameter 
limit, three more sightlines with a Lya absorber pass be¬ 
tween 200 and 300 kpc from a galaxy, which is between 
1.7 and 2 virial radii. In a very lenient approach, at best 
these might be classified as sampling the outermost halo 
of those galaxies. 

(5) We use the distribution of galaxies to define a 
filament at cz~3500 kms -1 . The filament has two 
branches with sizes 18°x7° (~15.5x6 Mpc) and 15°x7° 
(~13x6 Mpc), covering an area of 170 Mpc 2 . 

(6) All Lya lines with 1V(H I)>10 13 cm -2 originate 
within 2.1 Mpc of the filament axis. This includes de¬ 
tections that would be classified as “void” detections in 
an approach based solely on comparing their location to 
that of galaxies. 

(7) There is a strong correlation between the equivalent 
width of the Lya lines and filament impact parameter. 

(8) Using the simulations as a guide, we apply 
an ionization correction to each detected Lya line 
and find an average total hydrogen column density of 
loglV(H)=19.35. Combining with the area of the fila¬ 
ment this implies a total baryonic mass of 5.2x 10 13 M 0 . 
We also derive the combined mass of the galaxies inside 
the filament, which is about 1.4x10 13 Mq. 

(9) There is a strong suggestion that the Lya linewidth 
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correlates with filament impact parameter. 

(10) The four broad Lya components occur only in 
sightlines that pass close to (within 400 kpc) the axis of 
the filament. Further, multi-component absorbers pref¬ 
erentially occur within about 1 Mpc of the axis. Al¬ 
though the number of sightlines in our sample is small, 
this strongly suggest an increase in temperature and/or 
turbulence near the filament axis. 

(11) We note that Narayanan et al. (2010) previously 
found that the BLA toward Mrk 290 containing O VI pro¬ 
vides evidence for warm highly ionized gas in the filament 
with iV(H)/iV(H I)~1.8xl0 5 and a very large total col¬ 
umn density of 1.4x 10 19 cm -2 . This is consistent with 
the column density expected from the simulations. 

8.3. Comparison with simulations 

(12) Using the locations of galaxies in simulations we 
can find and define galaxy filaments that resemble the 
ones seen on the real sky. These filaments are then also 
seen in the H I column density distribution. 

(13) The predicted distribution of H I column den¬ 
sity as function of filament impact parameter is strongly 
dependent on the assumed intensity of the extragalactic 
background radiation field (EGB). Of the four versions of 
the Haardt & Madau (2001, 2012) EGB that we tested, 
we find that the best fit to the data occurs when using 
their standard 2001 version, in which the contribution to 
the ionizing flux at z=0 is about the same for quasars and 
galaxies. This implies that the HM12 model underesti¬ 
mates the ionizing flux by a factor four to five, in agree¬ 
ment with Kollmeier et al. (2014). Which is about the 
factor proposed by Khaire & Srianand (2015), but larger 
than the factor two proposed by Shull et al. (2015) to 
match the column density distribution of Lya absorbers. 

(14) Using the best matching EGB, the frac¬ 
tion of sightlines toward which we see Lya with 
IV(H I)> 10 13 cm -2 matches the simulations only at small 
filament impact parameter. For large filament impact 
parameters (>2.1 Mpc) the simulations predict a detec¬ 
tion rate of about 20-40% to find Lya absorbers with 
IV(H I)>10 13 cm~ 2 . The fact that we see 0 absorbers 
in 7 sightlines may be significant. Given the expected 
detection rate, probability of finding 0 detections in 7 
sightlines is 3-20%. This suggests that there may be a 
problem with the width of the filaments predicted by the 
simulations and/or the ionization background is stronger 
than HM01, though more data is needed. 

In summary, we have shown that the properties of the 
majority of Lya absorbers can be understood more easily 
with reference to the large scale structure of the Cosmic 
Web rather than to individual galaxy halos. By ana¬ 
lyzing the three-dimensional distribution of Lya forest 
absorbers it is possible to constrain simulations, as well 
as the extragalactic radiation field. More filaments need 
to be mapped to make our conclusions more statistically 
robust. 
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